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Abstract 

This  report  is  the  first  in  a  series  of  reports  giving  guidance  and  providing 
tools  to  practicing  engineers  on  the  proper  use  of  transient  seepage  analysis 
using  the  finite  element  method.  This  report  uses  innovative  analytical 
solutions  and  a  cross  section  of  a  generic  levee,  common  to  the  southeastern 
United  States,  to  provide  a  basic  introduction  to  properly  performing 
transient  seepage  analyses.  A  comparison  of  the  results  obtained  from  well- 
known  seepage  computer  programs  for  one-dimensional  analytical 
solutions  and  the  two-dimensional  generic  levee  cross  section  is  given. 

The  soil  property  data  and  geometry  data  needed  to  perform  a  transient 
seepage  analysis  are  discussed  in  detail.  Advice  on  when  a  transient  seepage 
solution  is  appropriate  is  given.  Hydrographs  that  mimic  typical  floods  are 
used  in  the  transient  analyses  for  the  generic  levee  cross  section.  Procedures 
for  comparing  the  results  of  different  transient  analyses,  as  well  as 
comparing  transient  analyses  to  steady  state  analyses,  are  developed. 
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Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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1  Introduction 

1.1  Description  of  need 

The  Risk  Management  Center  (RMC),  a  center  of  expertise  under  the 
U.S.  Army  Corps  of  Engineers  (USACE),  has  identified  a  research  need  to 
provide  more  guidance  on  performing  a  transient  seepage  analysis.  Transi¬ 
ent  seepage  analyses  have  become  an  ever-increasing  part  of  levee  analysis 
in  the  past  ten  years.  Many  commercial  computer  programs  are  available 
to  perform  these  types  of  analyses,  and  the  results  of  these  analyses  are 
becoming  much  more  common  in  engineering  reports.  The  USACE  has 
been  a  leader  in  developing  seepage  analysis  methods  and  design  guidance 
in  the  geotechnical  engineering  profession  for  the  past  8o  years,  and 
USACE  has  been  active  in  studying  the  mechanics  and  use  of  transient 
seepage  analysis  for  many  years  (Desai  1970;  Desai  1973;  Tracy  1973a  and 
1973b).  However,  the  implementation  of  transient  seepage  methods  by 
USACE  engineers  and  contractors  in  conventional  engineering  practice 
has  outpaced  the  development  of  guidance  documents  and  design 
recommendations. 

The  major  guidance  documents  used  for  seepage  analysis  of  levees  and 
dams  for  Corps  projects  are  Engineer  Manual  (EM)  1110-2-1901  (Seepage 
Analysis  and  Control  for  Dams)  (USACE  1993)  and  EM  1110-2-1913 
(Design  and  Construction  of  Levees)  (USACE  2000).  Although  these  EMs 
acknowledge  the  existence  of  transient  seepage  conditions,  they  do  not 
provide  guidance  regarding  the  use  of  transient  analyses  in  engineering 
practice. 

1.2  Objective  of  this  report 

The  objective  of  this  report  and  subsequent  work  remaining  on  this  project 
is  to  provide  information  to  help  guide  engineering  practitioners.  This 
report  is  the  first  submittal  as  part  of  a  larger  study  conducted  at  the 
U.S.  Army  Engineer  Research  and  Development  Center  (ERDC)  to  assess 
the  use  of  transient  seepage  analyses  in  Corps  of  Engineers  practice.  The 
near-term  questions  that  this  report  addresses  are  listed  below. 

1.  What  is  the  difference  between  a  transient  seepage  analysis  and  a  steady- 
state  seepage  analysis? 
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2.  What  additional  soil  parameters  are  required  for  conducting  a  transient 
seepage  analysis  as  compared  to  a  steady-state  seepage  analysis?  How  are 
these  parameters  obtained? 

3.  What  are  the  benefits  or  utility  of  conducting  transient  seepage  analyses? 
When  should  a  transient  analysis  be  used  in  lieu  of  a  steady-state  analysis? 

4.  What  computer  programs  can  be  used  for  transient  seepage  analysis?  How 
can  the  results  from  these  programs  be  verified? 

5.  What  factors  influence  the  results  of  transient  seepage  analysis? 

1.3  Approach  and  outline  of  tasks 

In  order  to  answer  these  questions,  the  following  procedure  was  adopted: 

1.  The  difference  between  steady-state  seepage  and  transient  seepage  is 
described. 

2.  The  equations  for  steady-state  seepage  and  transient  seepage  are 
developed  in  detail. 

3.  The  computer  programs  used  in  this  study  are  introduced. 

4.  Specific  issues  regarding  finite  element  analysis  of  transient  seepage  are 
discussed. 

5.  The  difference  in  the  soil  engineering  parameters  required  for  steady-state 
versus  transient  seepage  are  described  in  detail,  with  special  consideration 
given  to  their  implementation  in  the  computer  analyses. 

6.  Solutions  for  one-dimensional  (i-D)  transient  seepage  are  developed  for 
saturated  and  unsaturated  conditions.  The  closed  form  solutions  are 
compared  with  finite  element  solutions. 

7.  A  two-dimensional  (2-D)  cross  section  representing  a  “generic”  levee  is 
developed  and  analyzed  with  the  computer  programs. 

8.  The  sensitivity  of  the  finite  element  results  to  the  soil  volumetric  com¬ 
pressibility  and  other  factors  are  examined. 

General  recommendations  regarding  the  use  of  transient  seepage  analyses 

are  summarized  at  the  end  of  this  report. 
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2  Steady-state  versus  Transient  Seepage 
Analysis 

2.1  Differences  between  steady-state  seepage  and  transient 
seepage  anaiysis 

A  steady-state  condition  for  a  given  seepage  quantity  such  as  hydraulic 
head,  flow  rate,  or  a  given  soil  hydraulic  property  occurs  when  that  quantity 
is  not  changing  with  time.  Typically,  all  these  quantities  are  dependent  on 
each  other  so  that  if  one  is  varying,  the  others  are,  too.  However,  there  can 
be  instances  where  some  quantities  are  changing  while  others  are  not.  The 
given  quantity  is  in  a  transient  status  when  it  is  changing  with  time. 

Some  references  refer  to  steady-state  seepage  as  a  “saturated”  flow 
condition  and  transient  seepage  as  a  “partially  saturated”  or  “unsaturated” 
flow  condition.  While  some  elements  of  this  description  may  be  correct,  it  is 
not  a  generally  correct  concept.  Consider  the  case  of  steady-state  seepage 
stability  analysis  of  a  dam,  as  described  in  EM  1110-2-1902  (Slope  Stability) 
(USAGE  2003).  An  example  cross  section  of  a  dam  is  shown  in  Eigure  2.1. 
Eor  this  cross  section,  pore  pressures  are  determined  from  the  flow  net 
shown  on  the  figure,  which  was  drawn  for  steady-state  seepage  conditions. 
The  soil  is  saturated  below  the  phreatic  surface  and  partially  saturated  at 
some  distance  above  the  phreatic  surface.  Pore  pressures  can  also  be  cal¬ 
culated  using  the  finite  element  method.  In  order  to  determine  the  position 
of  the  phreatic  surface,  permeability  for  the  soil  for  both  saturated  and 
unsaturated  conditions  must  be  provided  to  the  finite  element  program. 
Elow  occurs  both  above  and  below  the  phreatic  surface.  This  is  an  example 
of  a  steady-state  flow  condition  where  the  soil  is  not  fully  saturated. 


Figure  2.1.  Cross  section  of  unconfined  steady-state  seepage  from  EM  1110-2-1901  (USAGE 

1993). 
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An  example  of  a  transient  flow  condition  where  the  soil  is  saturated  is  the 
falling  head  hydraulic  conductivity  test  (Figure  2.2).  This  test  is  used  to 
measure  the  hydraulic  conductivity  of  saturated  soils.  The  burette  attached 
to  the  top  of  the  saturated  sample  is  filled  with  water,  and  the  elapsed  time 
required  for  the  burette  to  drop  to  a  new  level  is  recorded.  The  flow 
through  the  sample  is  transient,  in  that  the  flow  decreases  with  time. 


Figure  2.2.  Schematic  of  a  faiiing  head  permeabiiity  test  iiiustrating 
transient  fiow  in  a  saturated  soii. 
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The  transient  seepage  case  of  most  interest  to  USAGE  is  when  the 
hydraulic  boundary  conditions  change  over  time.  When  a  steady-state 
seepage  analysis  is  performed  on  a  dam  or  levee,  it  is  assumed  that  the 
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flood  side  or  riverside  water  level  is  at  the  current  elevation  for  an  infinite 
amount  of  time.  This  is  shown  in  Figure  2.3a.  In  other  words,  time  is  not  a 
component  of  the  analysis  because  the  flow  is  constant.  Figure  2.3b  shows 
an  example  of  how  the  boundary  conditions  may  change  for  a  transient 
analysis.  At  an  elapsed  time  of  1  day,  there  may  be  no  water  on  the  levee. 
After  10  days  into  the  flood  event,  the  water  may  be  halfway  up  on  the 
levee.  The  water  level  may  continue  to  rise  to  a  higher  level  after  20  days, 
and  then  the  floodwater  may  have  receded  by  an  elapsed  time  of  30  days. 
This  is  the  type  of  hydraulic  loading  that  can  be  represented  by  a  transient 
seepage  analysis. 


Figure  2.3.  Steady-state  and  transient  boundary  conditions  on 
riverside  of  ievee. 


(a)  Steady  State  Analysis 


t  =  1  day 


t  =  20  days 


(b)  Transient  Analysis 
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2.2  Reasons  for  performing  a  transient  seepage  anaiysis 

The  mechanics  behind  transient  seepage  analyses  were  developed  more 
than  8o  years  ago.  However,  it  is  difficult  to  put  the  mechanics  into 
immediate  practical  use  owing  to  the  large  number  of  numerical  calcula¬ 
tions  required.  Only  after  robust  computer-based  numerical  modeling  tools 
were  developed,  such  as  finite  element  and  finite  difference  methods,  could 
transient  seepage  analyses  be  conducted  on  realistic  engineering  problems. 
In  recent  years,  transient  seepage  analysis  programs,  such  as  SEEP2D1, 
SEEP/W,  and  SLIDE,  allowed  these  analyses  to  be  conducted  on  personal 
computers  in  relatively  modest  execution  times.  However,  technology 
behind  performing  the  analysis  seems  to  have  outstripped  the  ability  to 
accurately  determine  reliable  input  parameters,  to  assess  the  validity  of  the 
results,  and  to  logically  incorporate  the  results  into  engineering  design. 

The  USAGE  has  been  conducting  steady-state  seepage  analyses  related  to 
levees  and  dams  for  more  than  70  years.  The  main  purposes  for  these 
analyses  have  included: 

1.  Determine  exit  vertical  gradients  to  calculate  factors  of  safety  against 
uplift. 

2.  Define  uplift  pressures  acting  on  top  blankets  of  low  permeability. 

3.  Calculate  flow  into  drains,  wells,  or  other  drainage  structures. 

4.  Determine  total  flow  beneath  structures. 

5.  Estimate  pore  pressures  for  effective  stress  slope  stability  analysis. 

Eor  the  analysis  cases  listed  above,  the  assumption  that  the  flow  has 
achieved  a  steady-state  condition  is  conservative  from  an  engineering 
perspective.  The  steady-state  analyses  will  normally  result  in  the  highest 
vertical  gradients,  uplift  pressures,  flows,  and  pore  pressures  that  the 
structures  should  experience.  There  has  been  a  growing  interest  in  sup¬ 
planting  these  steady-state  analyses  with  transient  analyses  for  various 
reasons,  with  some  being  justifiable  and  some  not. 

The  authors  have  seen  transient  analyses  included  in  engineering  reports 
dealing  with  routine  seepage  analyses  of  levees.  Although  these  transient 
analyses  were  not  requested  (nor  warranted),  it  appears  that  they  were 
conducted  because  (a)  the  computer  program  used  by  the  engineer  had  an 


1  Only  a  steady-state  version  of  SEEP2D  is  available  in  the  Groundwater  Modeling  System  (GMS)  suite  of 
programs.  A  transient  version  exists  and  has  been  used  internally  in  ERDC  projects.  As  part  of  this 
project,  a  transient  version  will  be  developed  for  GMS. 
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option  that  allowed  transient  analyses,  and  (b)  the  engineer  felt  that 
performing  transient  analyses  somehow  added  value  or  validity  to  their 
results.  In  most  of  these  cases,  the  transient  seepage  analyses  were 
performed  incorrectly  and  the  results  were  misleading. 

The  authors  have  also  seen  transient  analyses  being  used  to  explain  the 
results  of  piezometer  readings  during  transient  flooding  events.  This  can 
be  an  appropriate  use  of  transient  seepage  analyses  and,  under  some 
circumstances,  can  provide  reliable  results. 

Transient  analyses  can  often  be  successfully  used  to  estimate  the  time 
required  for  steady-state  seepage  conditions  to  be  achieved.  For  the  case  of 
determining  exit  gradients  for  the  factor  of  safety  against  uplift,  or  the 
heave  pressures  acting  on  the  base  of  a  top  stratum,  transient  analyses  can 
be  useful  to  estimate  the  development  of  the  uplift  forces  relative  to  the 
hydrograph  for  the  flood  event.  It  should  be  stressed  that  this  type  of 
analysis  only  provides  an  estimate.  It  may  not  produce  the  requisite 
accuracy  to  obtain  the  change  in  factor  of  safety  with  time. 

Transient  analyses  have  also  been  incorrectly  used  in  stability  analyses  to 
calculate  the  factor  of  safety  for  slope  stability  as  a  function  of  time  given  a 
flood  hydrograph.  In  this  type  of  analysis,  pore  pressures  from  transient 
seepage  analyses  have  been  paired  with  effective  stress  shear  strength 
parameters  to  determine  factors  of  safety.  This  type  of  analysis  assumes 
that  all  shear-induced  pore  pressures  have  been  dissipated.  For  the  case  of 
sand  levees,  this  assumption  may  be  valid.  For  the  case  of  levees 
constructed  of  fine-grained  soil,  this  assumption  may  not  be  correct,  and 
can  be  either  conservative  or  non-conservative.  Unfortunately,  this  type  of 
analysis  has  been  recommended  as  a  viable  analysis  method  in  a  variety  of 
published  sources  (Geo-Slope  2012;  Fredlund  et  al.  2011;  and  USER  2011). 

In  order  to  evaluate  the  different  programs  available  to  USAGE  personnel 
to  perform  transient  seepage  analysis,  and  to  examine  certain  cases  where 
a  transient  analysis  may  be  beneficial,  the  following  types  of  problems  are 
analyzed  in  this  report. 

1.  i-D  vertical  seepage  for  saturated  and  unsaturated  conditions. 

2.  2-D  steady-state  and  transient  seepage  for  a  hypothetical  levee  cross 
section. 

The  analysis  of  the  hypothetical  levee  cross  section  was  assessed  by 
determining  the  following  parameters: 
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1.  Total  head  values  at  specified  points  in  the  saturated  and  partially 
saturated  zones  as  a  function  of  time. 

2.  Exit  gradients  as  a  function  of  time  at  the  toe  of  the  levee. 

3.  Uplift  pressures  as  a  function  of  time  at  the  toe  of  the  levee. 

4.  Time  to  achieve  various  percentages  of  the  steady-state  conditions  for  the 
above  parameters. 

2.3  Computer  programs  used  in  study 

Three  computer  programs  were  used  as  part  of  this  study.  Two  are 
commercially  available  (SEEP/W  and  SLIDE),  and  one  is  under 
development  (SEEP2D).  SEEP/W  is  part  of  the  GeoStudio  2012 
(Version  8.12.2.7663)  suite  of  programs  produced  by  Geo-Slope 
International,  Ltd.  of  Calgary,  Alberta,  Canada.  It  is  a  stand-alone  finite 
element  seepage  program  that  can  do  steady-state  and  transient  seepage 
analysis.  Since  SEEP/W  is  an  unofficial  standard  for  almost  all  of  the 
Corps  of  Engineers,  it  is  used  for  all  the  computations  in  this  report  except 
where  comparisons  of  results  among  different  programs  are  made.  SLIDE 
(Version  6.024)  is  a  limit  equilibrium  slope  stability  program  that  is 
marketed  by  Rocscience,  Inc.,  of  Toronto,  Canada  (Rocscience  2010).  It 
also  contains  a  finite  element  seepage  model  that  can  do  steady-state  and 
transient  seepage  analyses. 

SEEP2D  was  derived  from  2-D  and  3-D  finite  element  computer  programs 
developed  for  the  USACE  Lower  Mississippi  Valley  Division  (LMVD)  to 
perform  steady-state  and  transient  seepage  analyses  (Tracy  1973a;  1973b). 
One  such  use  of  these  early  programs  was  a  2-D  and  3-D  seepage  analysis  of 
Lock  and  Dam  26  (Hall  et  al.  1975),  also  performed  for  LMVD.  As  the  state- 
of-the-art  continued  to  improve  and  because  the  overwhelming  majority  of 
seepage  analyses  by  USACE  was  2-D  steady  state,  an  improved  2-D,  steady- 
state  version  of  the  original  2-D  program  (Tracy  1973b)  was  developed 
(Tracy  1983).  This  program  became  known  as  SEEP2D.  A  seepage  package 
(Biedenharn  and  Tracy  1987)  was  later  developed  that  contained  a  mesh 
generation  program  (Tracy  1977b)  and  a  post-processor  program  (Tracy 
1977a)  for  plotting  results.  One  significant  improvement  in  SEEP2D  was  the 
addition  of  the  capability  to  generate  and  plot  a  flow  net.  In  1998,  SEEP2D 
was  incorporated  into  Version  2.1  of  the  Groundwater  Modeling  System 
(GMS  2013a;  2013b)  with  primary  work  done  by  Jones  (1999).  GMS, 
version  9.0.2,  is  currently  available  free  of  charge  to  USACE  personnel. 
Commercial  versions  are  available  from  Aquaveo  (2014). 
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3  Equations  for  Seepage 


The  governing  equations  for  seepage  flow  will  be  presented  with  the 
nomenclature  as  defined  in  Geo-Slope  (2012)  and  used  in  SEEP/W. 

3.1  Seepage  through  saturated  soils 

Equations  for  steady-state  and  transient  seepage  were  developed  for 
saturated  soils  and  presented  in  the  following  sections.  The  steady-state 
condition  is  the  simplest  case,  so  it  is  described  first. 


3.1.1  Steady-state  case 

Eor  the  case  of  in-plane  2-D  steady-state  flow  in  an  incompressible, 
homogeneous,  and  isotropic  soil  without  a  source  or  sink  such  as  a  well, 
the  governing  equation  is  the  well-known  Laplace’s  equation. 


dx^  dz^ 


(3.1) 


where: 

hf-  =  total  head  (L) 

X  =  X  coordinate  (L) 
z  =  z  coordinate  (L) 

When  the  soil  is  no  longer  homogeneous  and  a  source  or  sink  is  added,  the 
flow  equation  becomes 


A 

dx 


kOh. 


dx 


A 

dz 


kSh. 


dz 


+  Q*=0 


(3.2) 


where: 

k  =  hydraulic  conductivity  (L/T) 

Q*  =  source  (-1-)  or  sink  (-)  (i/T) 

When  the  soil  is  inhomogeneous  and  anisotropic  with  the  principal  axes  of 
flow  being  the  same  as  the  x-z  axes,  the  equation  of  flow  becomes 
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A 

dx 


d_ 

dz 


k  k 

ratio  X  o 

az 


+(?* 


0 


(3.3) 


where: 

kx  =  hydraulic  conductivity  in  the  first  principal  direction  (L/T) 
=  hydraulic  conductivity  in  the  second  principal  direction 
(L/T) 

kratio  ~  k^/kx 


Sometimes,  the  principal  axes  for  hydraulic  conductivity  do  not  coincide 
with  the  {x,z)  axes  (Figure  3.1)  but  are  offset  by  an  angle,  a.  In  this  final 
case,  the  steady-state  flow  equation  is 


2  7  *  2  \  ^  7  \  • 

cos  a +  sm  a)^  +  (l-A:,„,.Jsmacosa^ 


(l-Ario)si 


dK 


sm  a  cos  a  — ^  -f  ( sin  a  -f  k^„^^  cos  a 
dx 


2 


dz 


(3.4) 


+Q*=0 


Figure  3.1.  Principal  directions  of  hydraulic  conductivity  and  their  orientation  to  the 

(x,  z)  axes. 


X 
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3.1.2  Transient  case 

For  the  transient  case  and  where  the  soil-water  complex  is  now 
compressible,  the  equation  for  flow  becomes 


2  7  -2  \9K  1.1  \  •  dhf 

cos  a  +  /c^„,.„sm  a)^  +  (l-/c^„,.Jsmacosa^ 


.  ,  N  .  dh.  /  .  2  7  2  \9K 

l-^™77o)sinacosa^  +  (sm  a  +  /c^„,.„cos  a]^ 


(3.5) 


d9  dh 

dt  "  "  dt 


where: 

7vi/  =  density  of  water  (W/L3) 

rriy  =  slope  of  the  volumetric  strain  versus  effective  stress  curve  (see 
Figure  4.4)  (LVW) 

9  =  volumetric  moisture  content  (unitless) 


The  slope  of  the  volumetric  strain  versus  effective  stress  curve,  nty,  is 
treated  as  a  constant  for  transient  flow  through  a  saturated  soil. 

3.2  Transient  seepage  through  partially  saturated  soils 

Equation  3.5  is  the  governing  equation  for  transient  seepage  through 
partially  saturated  soils,  except  that  often  rriy  is  replaced  with  which  is 
the  negative  of  the  slope  of  the  volumetric  water  content  versus  pore  water 
pressure  curve.  In  the  case  of  partially  saturated  flow,  and  become 
functions  of  pore  water  pressure  and  therefore  must  be  calculated  from 
supplied  hydraulic  conductivity  and  volumetric  water  content  curves  as 
explained  in  Chapter  4. 
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4  Soil  Parameters  for  Seepage  Analyses 

The  soil  parameters  required  for  seepage  analyses  depend  on  the  relative 
complexity  of  the  analysis  conducted.  The  types  of  analyses,  listed  in  terms 
of  increasing  complexity,  are  shown  below: 

1.  Saturated  steady-state  confined  flow 

2.  Partially  saturated  steady-state  unconfined  flow 

3-  Saturated  transient  confined  flow 

4-  Partially  saturated  transient  unconfined  flow 

4.1  Saturated  steady-state  confined  flow 

Historically,  most  seepage  analyses  conducted  for  levees  have  been  for 
saturated  steady-state,  confined  flow  conditions.  In  these  types  of 
analyses,  the  top  flow  line  is  defined  by  the  geometry  of  the  problems  as 
opposed  to  being  controlled  by  the  soil  properties.  Shown  in  Figure  4.1 
below  is  an  example  of  this  type  of  problem.  The  levee  itself  is  assumed  to 
be  impervious.  This  is  consistent  with  the  Blanket  Theory  solutions  that 
have  been  used  by  the  USAGE  for  underseepage  calculations  over  the  past 
60  years.  The  top  flow  line  is  coincident  with  the  interface  between  the 
levee  and  the  foundation  soil. 


Figure  4.1.  Example  of  saturated  steady-state  confined  flow  problem  for  seepage  beneath 

levee. 


Pervious  Layer  (aquifer^ 
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For  this  type  of  problem,  the  main  soil  parameters  needed  are  the  satu¬ 
rated  hydraulic  conductivity  or  permeability  of  the  top  blanket  and  the 
pervious  layer.  No  other  soil  properties  are  required. 

4.2  Partially  saturated  steady-state  unconfined  flow 

A  more  complex  analysis  is  required  if  the  levee  is  pervious  relative  to  the 
top  blanket.  If  flow  is  occurring  through  the  levee,  then  the  top  flow  line  is  a 
function  of  the  soil  properties  of  the  levee,  and  the  flow  regime  is  “uncon¬ 
fined.”  This  type  of  flow  can  also  be  considered  “partially  saturated”  because 
the  soil,  for  some  distance  above  the  phreatic  surface  (line  of  zero  pore 
pressure),  is  partially  saturated,  and  the  soil  below  the  phreatic  surface  is 
considered  to  be  saturated.  A  schematic  of  this  type  of  seepage  problem  is 
shown  in  Figure  4.2.  For  this  category  of  problem,  the  saturated  perme¬ 
ability  values  for  the  top  blanket  and  pervious  layer  are  still  required,  but 
additionally,  a  hydraulic  conductivity  function  is  required  for  the  levee 
material.  The  hydraulic  conductivity  function  describes  the  value  of 
hydraulic  conductivity  as  a  function  of  the  matric  suction  in  the  soil  or  as  a 
function  of  the  volumetric  water  content  of  the  soil.  Matric  suction,  x(j,  is 
defined  as  the  difference  between  the  pore  air  pressure,  Ua,  and  the  pore 
water  pressure,  Uw  An  example  of  a  hydraulic  conductivity  function  is 
shown  in  Figure  4.3. 


Figure  4.2.  Example  of  unsaturated  steady-state  unconfined  flow  problem  for  seepage 

beneath  levee. 
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Figure  4.3.  Hydraulic  conductivity  functions. 


Matric  Suctiofi  ^f} 


It  is  difficult  to  accurately  determine  a  complete  hydraulic  conductivity 
function  for  a  soil.  It  is  relatively  common  to  conduct  laboratory  tests  to 
determine  the  saturated  hydraulic  conductivity  of  levee  fill  material.  The 
American  Society  for  Testing  and  Materials  (ASTM)  D5084  (Standard  Test 
Methods  for  Measurement  of  Hydraulic  Conductivity  of  Saturated  Porous 
Materials  Using  a  Flexible  Wall  Permeameter)  is  an  established  test 
procedure  and  can  be  conducted  by  most  geotechnical  testing  laboratories. 
It  is  much  more  difficult  to  determine  the  hydraulic  conductivity  for  clayey 
soils  at  elevated  values  of  matric  suction  and  for  partially  saturated  soils. 
There  is  a  testing  protocol  specified  by  ASTM  D7664  (Standard  Test  Meth¬ 
ods  for  Measurement  of  Hydraulic  Conductivity  of  Unsaturated  Soils),  but 
it  has  not  found  widespread  use  in  geotechnical  practice.  Few  commercial 
testing  laboratories  are  able  to  conduct  the  tests  described  in  the  ASTM 
standard.  In  many  cases,  it  is  possible  to  estimate  the  hydraulic  conductiv¬ 
ity  function  from  the  Soil  Moisture  Retention  Curve  (SMRC)  using  one  of 
an  assortment  of  published  procedures.  This  is  described  in  more  detail  in 
a  later  portion  of  this  report. 

The  required  accuracy  of  the  hydraulic  conductivity  function  depends 
upon  the  purpose  of  the  seepage  analysis.  If  the  analysis  is  being 
conducted  to  determine  the  vertical  exit  hydraulic  gradient  at  the  toe  of 
the  levee,  the  total  seepage  underneath  the  structure,  or  the  uplift 
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pressures  on  the  top  stratum,  then  an  approximate  function  for  the  levee 
fill  material  is  normally  adequate.  The  approximate  function  may  be  as 
simple  as  shown  as  Curve  a  in  Figure  4.3,  where  the  hydraulic  conductivity 
rapidly  decreases  as  the  pressure  head  goes  from  positive  to  negative.  The 
hydraulic  conductivity  for  all  values  of  negative  pressure  head  (positive 
matric  suction)  can  be  assumed  to  be  a  constant  factor  (0.001  for  example) 
multiplied  by  the  saturated  hydraulic  conductivity  or  permeability  of  the 
levee  fill  material.  Alternatively,  generic  hydraulic  conductivity  functions 
are  available  via  pull-down  menus  for  various  soil  types  in  the  seepage 
analysis  software.  Curve  b  shown  in  Figure  4.3  is  the  function  provided  by 
SLIDE  for  a  soil  type  of  “clay.” 

If  the  purpose  of  the  seepage  analysis  is  to  determine  pore  pressure  values 
in  the  levee  fill  for  use  in  an  effective  stress  steady-state  stability  analysis, 
such  as  that  described  in  EM  1110-2-1902,  then  an  accurate  hydraulic 
conductivity  function  is  desired.  The  hydraulic  conductivity  function 
influences  the  position  of  the  phreatic  surface,  which,  in  turn,  influences 
the  values  of  pore  pressure.  In  simple  terms,  the  choice  of  the  hydraulic 
conductivity  function  (e.g..  Curve  a  versus  Curve  b  in  Eigure  4.3)  can 
influence  the  calculated  stability  factor  of  safety. 

There  are  many  different  methods  to  determine  more  refined  hydraulic 
conductivity  functions  than  those  shown  in  Eigure  4.3  that  are 
incorporated  into  commercial  seepage  software.  Most  of  the  methods  used 
in  the  software  were  developed  from  research  linking  the  water  content 
function  with  the  hydraulic  conductivity  function,  and  several  of  these 
methods  are  listed  in  Table  4.1.  In  SEEP/W,  the  water  content  function 
must  first  be  specified  to  obtain  the  hydraulic  conductivity  function.  In 
SLIDE,  the  hydraulic  conductivity  function  can  be  determined  for  steady- 
state  seepage  problems  without  specifically  defining  the  water  content 
function.  All  methods  require,  as  a  minimum,  the  value  of  the  saturated 
hydraulic  conductivity. 


Table  4.1.  Methods  to  estimate  the  hydraulic  conductivity  function  in  SLIDE,  SEEP2D,  and 

SEEP/W. 


Method 

SEEP2D1 

SEEP/W 

SLIDE 

Brooks  and  Corey  (1964) 

Fredlund  and  Xing  (1994) 

Gardner  (1958) 

van  Genuchten  (1980) 

1  other  methods  to  estimate  the  hydraulic  conductivity  function  are  planned  to  be  added  to  SEEP2D. 
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Each  of  the  different  methods  has  certain  advantages  and  disadvantages 
for  specific  soil  types.  It  is  beyond  the  scope  of  this  initial  report  to  make 
recommendations  regarding  the  choice  of  the  specific  method  to  be  used. 

A  subsequent  report  will  be  submitted  that  addresses  the  use  of  the  differ¬ 
ent  methods. 

4.3  Saturated  transient  confined  fiow 

A  saturated  transient  confined  flow  seepage  problem  can  be  performed  on 
the  cross  section  shown  in  Figure  4.1.  If  the  water  level  on  the  riverside  of 
the  levee  varies  with  time,  then  the  flow  beneath  the  levee  would  also  vary 
with  time,  and  this  would  constitute  a  transient  seepage  condition.  Along 
with  specifying  the  change  in  hydraulic  boundary  conditions  with  time,  the 
changes  in  the  soils’  hydraulic  properties  with  pore  water  pressure  are 
required  in  the  analysis. 

Because  the  top  stratum  and  the  pervious  layer  are  both  saturated,  and 
will  remain  saturated  during  the  flow  event,  a  hydraulic  conductivity  func¬ 
tion  is  not  necessary.  However,  a  value  of  the  coefficient  of  volumetric 
compressibility,  mv,  is  required.  This  parameter  is  the  slope  of  a  tangent 
drawn  to  a  point  on  a  conventional  consolidation  curve  for  an  arithmetic 
stress  axis  as  shown  in  Figure  4.4.  The  magnitude  of  the  coefficient  of 
volumetric  compressibility  has  a  very  significant  effect  on  the  results  of 
transient  seepage  analyses.  This  will  be  discussed  later  in  the  report. 
Approximate  values  of  are  given  in  Table  4.2. 


Figure  4.4.  Coefficient  of  voiumetric  compressibiiity  shown  for 
conventionai  consoiidation  test  data. 


Effective  Stress 
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Table  4.2.  Range  of  nriv  values  for  various  soils  (after  Domenico  and  Mifflin  1965). 


Soil  Type 

mv{VPa-^) 

/77i.(psU) 

Plastic  clay 

2.1x10-3  to  2.6x10-4 

1x10-4  to  1.25x10-5 

Stiff  clay 

2.6x10-4  to  1.3x10-4 

1.25x10-5  to  6.25x10-6 

Medium  hard  clay 

1.3x10-4  to  6.9x10-5 

6.25x10-6  to  3.3x10-6 

Loose  sand 

lx  10-4  to  5.2x10-5 

5x10-6  to  2.5x10-6 

Dense  sand 

2.1x10-5  to  1.3x10-5 

1x10-6  to  6.25x10-7 

Dense  sandy  gravel 

1x10-5  to  5.2x10-6 

5x10-7  to  2.5x10-7 

Jointed  rock 

6.9x10-6  to  3.3x10-7 

3.3x10-7  to  1.6x10-8 

Sound  rock 

>3.3x10-7 

>1.6x10-8 

Water 

4.4x10-7 

2.1x10-8 

SEEP/W  has  a  special  “saturated  only”  soil  model  that  can  be  used  for  this 
seepage  case.  In  addition  to  the  value  of  saturated  hydraulic  conductivity 
and  compressibility  values,  SEEP/W  also  requires  the  user  to  input  the 
value  of  the  saturated  volumetric  water  content,  0^,  for  the  soil  layers.  Eor 
all  transient  seepage  analysis  cases,  SLIDE  requires  a  complete  set  of  soil 
properties  to  be  entered  (e.g.,  water  content  function  and  hydraulic  con¬ 
ductivity  function),  and  it  does  not  have  a  special  “saturated  only”  soil 
model.  Transient  SEEP2D  has  not  yet  been  finalized,  but  it  will  be  similar 
to  SEEP/W  and  SLIDE. 

4.4  Partially  saturated  transient  unconfined  flow 

The  final  seepage  analysis  case  is  similar  to  that  shown  in  Eigure  4.2,  except 
that  the  riverside  water  level  changes  as  a  function  of  time.  Elow  occurs 
through  the  levee  fill  material  as  it  becomes  more  and  more  saturated  over 
time.  After  enough  time  has  passed,  a  steady-state  flow  condition  will  be 
achieved,  and  the  phreatic  surface,  as  shown  in  Eigure  4.2,  will  have  been 
developed.  Eor  this  seepage  case,  a  complete  set  of  property  data  is  required 
for  the  levee  fill  material.  Since  the  top  stratum  and  the  pervious  layer 
would  remain  saturated,  the  “saturated  only”  soil  model  can  be  used  for 
these  materials  in  SEEP/W.  If  SLIDE  is  used,  a  complete  set  of  property 
data  is  required  for  all  soil  layers. 

A  complete  set  of  property  data  for  transient  seepage  cases  includes  (1)  the 
saturated  permeability,  (2)  the  water  content  function,  (3)  the  hydraulic 
conductivity  function,  and  (4)  the  coefficient  of  volumetric  compressibility. 
The  saturated  permeability  is  a  part  of  all  seepage  analyses  and  is  covered  in 
detail  in  EM  1110-2-1901  (USAGE  1993)  and  other  USAGE  manuals  (USAGE 
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1956;  USAGE  1986;  USAGE  1998;  USAGE  2000;  USAGE  2005).  Details 
regarding  the  water  content  function  and  hydraulic  conductivity  function 
are  further  discussed  in  this  section. 

4.5  Water  content  function 

The  water  content  function  characterizes  the  decrease  in  moisture  content 
as  the  matric  suction  increases  on  a  soil  sample.  As  the  suction  increases, 
water  is  pulled  from  the  soil,  and  the  water  content  decreases.  In  unsatu¬ 
rated  soil  mechanics,  the  volumetric  water  content  (0)  is  most  often  used. 
This  is  defined  as  the  volume  of  water  divided  by  the  total  volume  of  the  soil 
sample.  Eor  a  saturated  soil,  the  volumetric  water  content  (0s)  is  numeri¬ 
cally  equal  to  the  porosity,  n.  It  is  important  not  to  confuse  the  volumetric 
water  content  used  in  transient  seepage  analyses  with  the  gravimetric  water 
content  (w)  that  is  normally  used  in  geotechnical  engineering  practice.  The 
water  content-  matric  suction  interrelationship  is  demonstrated 
schematically  in  Eigure  4.5.  If  the  level  of  the  water  in  the  burette  is  at  the 
same  elevation  as  the  top  of  the  soil  sample,  and  equilibrium  has  been 
achieved,  the  matric  suction  is  equal  to  zero,  and  the  pore  water  pressure  is 
equal  to  zero.  The  volumetric  water  content  of  the  soil  sample  would  be 
equal  to  Og  at  this  point.  If  the  elevation  of  the  burette  is  lowered,  the  matric 
suction  is  increased,  the  pore  pressure  is  decreased,  and  water  flows  from 
the  soil  sample  into  the  burette.  The  soil  sample  will  reach  a  new 
equilibrium  water  content  as  a  function  of  the  matric  suction  applied.  If  the 
volumetric  water  content  is  measured  for  a  range  of  matric  suction,  the 
water  content function  can  be  plotted.  An  example  water  content  function 
is  shown  in  Eigure  4.6.  The  computer  programs  investigated  as  part  of  this 
study  allow  the  water  content  functions  to  be  directly  entered  as  parameter 
values.  It  is  important  to  note  that  the  water  content  function  depends  on  if 
the  soil  is  drying  (desorption)  or  wetting  (adsorption).  Normally,  only  one 
water  content  function  can  be  input  for  each  soil,  and  the  engineer  needs  to 
make  a  choice  on  which  curve  (drying  versus  wetting)  to  use,  or  if  an 
average  water  content  function  should  be  used. 
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Figure  4.5.  Schematic  of  test  to  increase  matric  suction  and  decrease  soii 

moisture  content. 


Figure  4.6.  Exampie  water  content  function  (Frediund  et  ai.  1994). 


Owing  to  the  multidisciplinary  efforts  involved  in  research  into  transient 
seepage  analysis,  there  can  be  some  confusion  regarding  the  nomenclature 
used  for  the  various  soil  parameters.  The  water  content  function  is  also 
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called  the  soil  moisture  retention  eurve  (SMRC),  soil  water  retention 
eurve  (SWRC),  water  retention  eharaeteristie  eurve  (WRCC),  storage 
funetion,  or  suetion  head  eurve.  When  the  curve  is  plotted,  one  axis  is 
normally  the  volumetric  water  content.  The  other  axis  may  be  either  the 
matric  suction  (xp)  or  the  pore  pressure  (u),  and  this  axis  may  be  plotted  as 
an  arithmetic  or  logarithmic  scale.  Shown  in  Figure  4.7  is  an  alternative 
method  of  plotting  the  water  content  function.  This  type  of  plot  clearly 
shows  how  the  coefficient  of  volume  compressibility  relates  to  the  water 
content  function.  Figure  4.8  shows  typical  water  content  functions  for 
sand,  silt,  and  clay. 

The  test  apparatus  shown  in  Figure  4.5  is  similar  to  a  Tempe  cell  test. 
There  is  a  limit  to  the  value  of  matric  suction  that  can  be  applied  using  an 
apparatus  of  this  sort,  but  other  test  apparatuses  are  available  to  measure 
the  water  content  function  for  high  values  of  matric  suction.  There  are 
many  different  tests  that  have  been  proposed  to  measure  the  water  content 
function  of  soils.  Many  of  these  tests  have  been  standardized  by  ASTM  and 
these  are  listed  in  Table  4.3.  This  is  currently  a  popular  research  area,  and 
new  test  methods  and  apparatuses  are  being  introduced  frequently. 

Figure  4.7.  Water  content  function  showing  the  coefficient  of  voiume 
compressibiiity  (after  Geo-studios  2013). 
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Figure  4.8.  Water  content  function  for  different  soii  types  (after  Geo¬ 
studios  2013). 
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Tabie  4.3.  Tests  for  transient  seepage  parameters  standardized  by  ASTM. 


Std.  No. 

Titie 

D2325 

Test  Method  for  Capillary-Moisture  Relationships  for  Coarse-  and  Medium- 
Textured  Soils  by  Porous-Plate  Apparatus 

D3152 

Test  Method  for  Capillary-Moisture  Relationships  for  Fine-Textured  Soils  by 
Pressure-Membrane  Apparatus 

D5298 

Standard  Test  Method  for  Measurement  of  Soil  Potential  (Suction)  Using  Filter 
Paper 

D6527 

Test  Method  for  Determining  Unsaturated  and  Saturated  Hydraulic  Conductivity 
in  Porous  Media  by  Steady-state  Centrifugation 

D6836 

Test  Methods  for  Determination  of  the  Soil  Water  Characteristic  Curve  for 
Desorption  Using  Hanging  Column,  Pressure  Extractor,  Chilled  Mirror 

Hygrometer,  or  Centrifuge 

D7664 

Standard  Test  Methods  for  Measurement  of  Hydraulic  Conductivity  of 

Unsaturated  Soils 
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4.6  Estimation  of  water  content  function 

The  computer  programs  investigated  as  part  of  this  research  (SEEP2D, 
SEEP/W,  and  SLIDE)  can  accept  the  water  content  functions  as  discrete 
coordinates.  Methods  are  available  to  estimate  the  curves  based  on  grain 
size  information,  soil  type,  or  previously  published  data.  SEEP/W  has 
“generic”  water  content  functions  for  six  different  soil  types.  Alternatively, 
the  water  content  function  can  be  estimated  from  gradation  data  using  the 
method  described  by  Aubertin  et  al.  (2003).  If  the  water  content  function 
is  input  as  discrete  coordinates,  then  the  Eredlund  and  Xing  (1994)  or  van 
Genuchten  (1980)  methods  can  be  used  to  fit  curves  to  the  data,  and  the 
resulting  curve-fit  parameters  can  be  used  for  prediction  of  the  hydraulic 
conductivity  function. 

SLIDE  allows  you  to  select  from  four  soil  types  (sand,  silt,  clay,  and  loam) 
or  to  use  one  of  the  four  methods  listed  in  Table  4.1  to  estimate  the  water 
content  function.  Eor  each  of  the  four  methods,  the  user  can  select  the 
appropriate  scaling  or  curve-fit  parameters  based  on  published  results  for 
specific  soils,  or  individual  parameters  can  be  entered. 

4.7  Estimation  of  hydrauiic  conductivity  function 

The  hydraulic  conductivity  functions  can  be  entered  as  discrete  coordinates 
if  data  are  available,  but  this  is  currently  rare  in  geotechnical  engineering 
practice.  The  hydraulic  conductivity  function  is  most  often  estimated  from 
the  water  content  function  using  one  of  the  methods  listed  in  Table  4.1.  It  is 
beyond  the  scope  of  this  initial  report  to  make  recommendations  of  the 
method  to  use,  but  the  most  common  methods  in  engineering  practice  seem 
to  be  those  described  in  van  Genuchten  (1980)  and  Eredlund  and  Xing 
(1994)- 

As  new  methods  are  developed,  they  can  be  incorporated  into  transient 
seepage  analysis  by  separately  calculating  the  water  content  function  and 
the  hydraulic  conductivity  function  in  a  spreadsheet,  and  then  inputting 
those  data  directly  into  the  computer  program.  As  an  example.  Sleep 
(2011)  proposed  a  method  to  estimate  the  water  content  function  based  on 
the  saturated  hydraulic  conductivity  of  the  soil.  This  method  can  be  used 
to  develop  a  table  of  the  volumetric  water  content  and  matric  suction 
values  that  can  by  input  into  the  computer  program.  The  hydraulic 
conductivity  function  can  then  be  calculated  from  these  data  using  the 
method  of  van  Genuchten  (1980)  or  Eredlund  and  Xing  (1994). 
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5  Basic  Concepts  of  the  Finite  Eiement 
Method  for  Seepage 

5.1  Finite  eiement  mesh 

The  programs  used  for  this  study  have  the  ability  to  automatically  generate 
the  finite  element  mesh.  SEEP2D  currently  uses  3-noded  triangular 
elements  with  4-noded  quadrilateral  element  capability  to  be  added  in  the 
next  version.  Both  SLIDE  and  SEEP/W  allow  the  use  of  3-noded  and 
6-noded  triangular  elements  and  4-noded  and  8-noded  quadrilateral 
elements.  SEEP/W  allows  triangular  and  quadrilateral  elements  to  be 
mixed  in  a  single  mesh.  As  a  general  rule,  the  relative  element  size  in  the 
finite  element  mesh  should  be  similar  to  the  relative  size  of  squares  in  a 
flow  net.  Smaller  elements  should  be  used  in  areas  of  high  gradients  and 
also  in  areas  of  particular  engineering  interest.  Additional  guidance 
regarding  the  construction  of  the  finite  element  mesh  for  seepage  analyses 
can  be  found  in  Duncan  et  al.  (2011). 

5.2  Convergence 

The  concept  of  convergence  of  the  finite  element  solution  will  now  be 
discussed.  A  successfully  converged  solution  is  critical  to  an  accurate 
simulation. 

5.2.1  Description  and  program  impiementation 

Whenever  there  is  unsaturated  flow,  the  computer  programs  must  iterate 
to  arrive  at  a  solution.  This  can  be  explained  in  various  ways.  Eirst,  from  a 
physical  point  of  view  and  for  the  steady-state  case,  the  positions  of  the 
phreatic  surface  and  exit  point  where  the  phreatic  surface  intersects  the 
exit  face  are  not  initially  known.  In  SEEP2D,  an  initial  “guess”  is  made, 
and  the  linear  system  of  equations  is  formulated  and  solved  for  the  total 
head  at  each  node  point.  Eor  all  elements  in  the  unsaturated  zone,  their 
respective  hydraulic  conductivities  will  be  changed  to  values  determined 
by  the  hydraulic  conductivity  function.  The  linear  system  of  equations  is 
then  reformulated,  and  total  head  at  each  node  is  recomputed.  The  com¬ 
puted  total  heads  at  the  nodes  will  almost  always  be  different  from  the 
initial  guess.  When  the  system  of  linear  equations  is  formulated  a  second 
time,  its  coefficients  for  nodes  in  the  unsaturated  zone  will  be  different 
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from  the  first  solution,  so  the  computed  total  head  at  each  node  for  this 
second  system  of  equations  is  different  from  the  first  solution.  This 
process  is  repeated  until  the  amount  of  change  of  total  head  at  all  the 
nodes  remains  very  close  when  comparing  the  last  two  solutions.  When 
this  process  has  achieved  an  acceptable  level  of  accuracy,  the  solution  is 
said  to  have  converged.  Mathematically,  Equation  3.3  governing  steady- 
state  flow  requiring  this  iterative  procedure  for  solution  is  said  to  be 
nonlinear. 

For  the  transient  case,  the  position  of  the  phreatic  surface  and  exit  point  is 
known  from  the  initial  conditions  of  the  problem  where  a  total  head  at  all 
the  node  points  is  provided.  Equation  3.5  representing  transient  flow  is 
also  nonlinear,  so  an  iteration  sequence  is  required  as  in  the  steady-state 
case  for  each  time  step  of  the  solution  process.  The  only  difference  in  the 
transient  case  is  the  linear  system  of  equations  resulting  from  Equation  3.5 
is  different. 

Fewer  details  are  available  about  the  convergence  and  solution  schemes 
used  by  the  commercial  software  packages.  However,  there  are  features  in 
the  programs  that  allow  the  user  to  determine  when  convergence  has  been 
achieved.  In  SEEP/W,  the  Draw-Graph  function  can  be  used  to  show  a 
plot  of  the  number  of  unconverged  pressure-head  nodes  for  each  time 
step,  or  the  results  can  be  shown  in  table  form.  In  addition,  convergence 
information  is  displayed  on  the  computer  screen  as  the  program  calculates 
the  solution.  SLIDE  also  displays  convergence  information  during 
program  execution,  and  it  can  generate  a  similar  convergence  plot  after 
execution  using  the  program  module  Slideinterpret.  During  execution, 
SEEP2D  does  a  simple  printout  of  important  information  such  as  the 
absolute  value  of  the  maximum  difference  among  all  the  nodes  between 
the  last  two  iterations  for  a  given  nonlinear  iteration. 

5.2.2  Procedures  to  aid  convergence 

Obtaining  convergence  with  the  computer  programs  was  a  problem  area  in 
this  study.  The  goal  is  to  have  no  non-converging  nodes.  It  is  not  possible 
to  provide  guidance  for  overcoming  convergence  problems  that  always 
works,  but  here  are  some  general  rules  that  seem  to  improve  convergence. 

•  Increasing  the  number  of  nodes  and  elements. 

•  Decreasing  the  size  of  the  time  step. 
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•  Using  lower  order  elements  (3-noded  triangles  and  4-noded  quadri¬ 
laterals)  as  opposed  to  higher  order  elements  (6-noded  triangles  and 
8-noded  quadrilaterals). 

•  Using  “smooth”  water  content  and  hydraulic  conductivity  functions. 
Increasing  the  number  of  discrete  points  representing  these  curves  can 
sometimes  help. 

•  Using  “simpler”  water  content  and  hydraulic  conductivity  functions. 
Using  curves  with  less  severe  curvature  may  degrade  the  quality  of 
solution. 

•  Using  a  less  stringent  convergence  criterion.  However,  this  again  may 
degrade  the  quality  of  solution.  In  SEEP2D,  a  maximum  allowable  dif¬ 
ference  in  total  head  among  all  the  nodes  between  the  final  two  time 
steps  is  used  to  determine  convergence.  A  typical  value  of  0.001  with 
some  nodes  not  converging  may  well  be  worse  than  using  a  value  of 
0.05  with  all  nodes  converging. 

Increasing  the  number  of  elements  and  decreasing  the  size  of  the  time 
steps  also  increases  the  execution  time  required  to  obtain  a  solution.  In 
some  cases,  the  increase  in  execution  time  can  be  significant.  Some  of  the 
runs  conducted  during  this  study  took  greater  than  10  hr. 
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6  1-D  Transient  Seepage  -  Green-Ampt 

Problem 

In  order  to  become  acquainted  with  the  operation  of  transient  seepage 
computer  programs,  two  simple  i-D  problems  are  presented.  These 
problems  can  be  used  by  the  analyst  to  become  familiar  with  the  basic 
elements  of  conducting  a  transient  analysis  (i.e.,  meshing,  entering  soil 
parameters,  setting  boundary  conditions,  etc.)  Analytic  solutions  for  these 
problems  are  provided  in  detail  in  Appendices  A  and  B. 

The  Green-Ampt  infiltration  problem  (Green  and  Ampt  1911;  Smith  2010) 
is  typically  i-D  transient  unsaturated  flow  in  a  vertical  column  of  soil  of 
length  =  L  that  is  almost  dry  with  only  residual  moisture  content. 
Infiltration  is  applied  at  the  top  of  the  column  with  the  total  head  =  L.z  - 
0  at  the  bottom  of  the  column  and  is  positive  upward.  A  schematic  of  the 
geometry  of  this  problem  is  shown  in  Figure  6.1.  The  original  Green-Ampt 
(1911)  solution  assumed  a  wetting  front  of  completely  saturated  soil 
proceeding  downward  as  the  soil  goes  from  almost  dry  to  completely 
saturated.  However,  the  actual  solution  is  smooth. 


Figure  6.1.  Column  of  soil. 
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6.1  Unsaturated  flow 

While  the  water  advances  to  the  bottom  of  the  column,  the  flow  is  not 
constant  over  time  since  the  soil  is  unsaturated.  This  is  often  called 
“unsaturated  flow.” 

6.1.1  Initial  conditions 

The  almost  dry  soil  at  t  =  0  is  modeled  by 

ht=Pr+z  (6.1) 


where: 


Pr  =  residual  pressure  head  (L)  that  has  a  large  absolute  value  but 
negative. 

6.1.2  Boundary  conditions 

The  boundary  condition  at  the  bottom  (z  =  0)  is  /i^  =  Pr,  and  at  the  the  top 
(z  =  L)  is  hf-  —  L. 

6.1.3  Hydraulic  conductivity  function 

The  hydraulic  conductivity  function  can  be  represented  by 

k  =  KK  (6.2) 


where: 

kg  =  saturated  hydraulic  conductivity  (L/T) 
ky  =  relative  hydraulic  conductivity  (unitless) 

To  derive  the  analytic  solution,  the  simpler  Gardner’s  equation  (Gardner 
1958)  for  ky  is  used.  It  is 


K  =  e‘‘P,p<0;k^  =  l,p>0 


(6.3) 


where: 


p  =  pressure  head  (L) 
a  =  positive  parameter  (i/L) 
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6.1.4  Moisture  content  function 

First,  the  effective  saturation  is  defined  as 

C  _ 


(6.4) 


where: 

9r  =  residual  moisture  content  (unitless) 
9s  =  saturated  moisture  content  (unitless) 

Again  for  simplicity, 


Therefore,  the  soil-moisture  function  is 

0  =  + 


(6.5) 


(6.6) 


6.1.5  Finite  eiement  anaiysis  soiution 

The  analytic  solution  given  in  Appendix  A  was  used  to  test  the  computer 
programs,  SEEP/W,  SLIDE,  and  SEEP2D.  Eigure  6.2  shows  the  3-noded 
triangular  mesh  used  by  SEEP2D.  The  SLIDE  analysis  used  6-  noded 
triangular  elements,  and  the  SEEP/W  analysis  used  4-noded  quadrilateral 
elements.  The  values  of  the  input  parameters  used  in  the  analysis  are  given 
in  Table  6.1. 

Table  6.2  gives  the  results  after  a  period  of  1  day  and  5  days,  and 
Eigures  6.3  and  6.4  show  plots  of  the  computed  results  and  analytic 
solution.  Considerable  effort  was  needed  to  get  the  three  different 
programs  to  match  results.  These  problems  were  often  related  to  the  time 
steps  used,  how  convergence  was  implemented,  and  what  type  of  finite 
element  was  used.  Once  these  problems  were  solved,  the  results  of  all 
three  programs  matched  closely  with  the  analytic  solution. 
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Figure  6.2.  Finite  eiement  mesh  with  boundary  conditions. 


Tabie  6.1.  Parameters  and  vaiues  used  in  the  unsaturated  fiow  exampie. 


L  =  50  ft 

Height  of  soil  column 

II 

CO 

o 

Residual  pressure  head  or  negative  the  matric  suction 

=  0.15 

Residual  water  content 

B,  =  0.45 

Saturated  water  content 

ks  =  0.1  ft/day 

Saturated  hydraulic  conductivity 

Az  =  2.0  ft 

Height  of  elements  in  z  direction 

Ax  =  2.0  ft 

Width  of  elements  in  x  direction 

At  =  0.001  days 

Time  step  (may  be  constant  or  variable) 

a  =  0.1  ft-i 

Constant  in  Gardner  (1958)  equation 

=  10-5  ftVIb 

Slope  of  the  volumetric  strain  versus  effective  stress  curve 
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Table  6.2.  Total  head  at  different  elevations  and  after  different  times  for  the  analytic  solution,  SEEP2D, 
SEEP/W,  and  SLIDE  for  the  unsaturated  flow  problem. 


Total  Head  (ft) 

Time  =  1  day 

Time  =  5  days 

z{ft) 

Analytic 

SEEP2D 

SEEP/W 

SLIDE 

Analytic 

SEEP2D 

SEEP/W 

SLIDE 

48 

41.24 

40.60 

41.01 

41.41 

45.87 

45.83 

45.84 

45.97 

46 

29.38 

28.87 

29.01 

29.83 

41.01 

40.90 

40.95 

41.15 

44 

18.16 

18.83 

18.14 

18.30 

35.43 

35.25 

35.34 

35.46 

42 

12.54 

13.16 

12.66 

12.43 

29.23 

29.00 

29.10 

28.93 

40 

10.03 

10.21 

10.07 

10.02 

22.65 

22.47 

22.52 

21.78 

38 

8.00 

8.03 

8.01 

8.00 

16.19 

16.16 

16.08 

14.73 

36 

6.00 

6.00 

6.00 

6.00 

10.51 

10.66 

10.44 

8.90 

34 

4.00 

4.00 

4.00 

4.00 

6.06 

6.35 

6.06 

4.92 

32 

2.00 

2.00 

2.00 

2.00 

2.79 

3.02 

2.81 

2.23 

30 

0.00 

0.00 

0.00 

0.00 

0.26 

0.41 

0.28 

0.05 

28 

-2.00 

-2.00 

-2.00 

-2.00 

-1.92 

-1.86 

-1.91 

-1.99 

26 

-4.00 

-4.00 

-4.00 

-4.00 

-3.98 

-3.95 

-3.97 

-4.00 

24 

-6.00 

-6.00 

-6.00 

-6.00 

-6.00 

-5.98 

-5.99 

-6.00 

22 

-8.00 

-8.00 

-8.00 

-8.00 

-8.00 

-8.00 

-7.99 

-8.00 

20 

-10.00 

-10.00 

-10.00 

-10.00 

-10.00 

-10.00 

-10.00 

-10.00 

6.2  Saturated  flow 

6.2.1  Description  of  the  probiem 

For  the  second  example,  the  column  of  soil  is  initially  fully  saturated,  having 
a  total  head  of  =  Hq  -  L  until  a  total  head  of  -  H^>  Hq  is  applied  at 
the  top  (z  =  L)  and  remains  at  this  level  for  an  indefinite  period  of  time.  The 
total  head  at  the  bottom  (z=  0)  remains  hf-  =  Hq  .  Total  head  then  increases 
down  the  column  until  the  steady-state  condition  is  achieved. 

The  initial  condition  is  that  no  flow  is  occurring  at  t  =  0  since  a  constant 
value  of  ht  =  Hq  exists  throughout  the  column.  The  boundary  condition  at 
z  =  0  is  /if  =  Hq,  and  at  z  =  L,  it  is 
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Figure  6.3.  Total  head  (ft)  for  different  z  values  (ft)  for  the  analytic, 
SEEP/W,  SLIDE,  and  SEEP2D  solutions  at  Time  =  1  day  for  the 
unsaturated  flow  problem. 


Figure  6.4.  Total  head  (ft)  for  different  z  values  (ft)  for  the  analytic, 
SEEP/W,  SLIDE,  and  SEEP2D  solutions  at  Time  =  5  days  for  the 
unsaturated  flow  problem. 
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6.2.2  Finite  eiement  anaiysis  soiution 

The  analytic  solution  given  in  Appendix  B  for  this  saturated  transient  flow 
problem  was  used  to  assess  the  computer  programs,  SEEP/W,  SLIDE,  and 
SEEP2D,  for  the  saturated  transient  solution.  The  same  finite  element 
meshes  were  used  for  these  analyses  as  used  for  the  unsaturated  analyses. 
The  values  of  the  input  parameters  used  in  the  analysis  are  given  in 
Table  6.3. 


Table  6.3.  Parameters  and  values  used  In  the  saturated  flow  example. 


Parameter  and  Value 

Description 

L  =  50  ft 

Height  of  soil  column 

Ho  =  50  ft 

Initial  total  head 

II 

CJl 

Final  total  head  at  top  of  column 

Or  =  0.15 

Residual  water  content 

Os  =  0.45 

Saturated  water  content 

ks  =  0.1  ft/day 

Saturated  hydraulic  conductivity 

Az  =  2.0  ft 

Height  of  elements  in  z  direction 

Ax  =  2.0  ft 

Width  of  elements  in  x  direction 

At  =  0.001  days 

Time  step  (may  be  constant  or  variable) 

=  10-5  ftVIb 

Slope  of  the  volumetric  strain  versus  effective  stress  curve 

Table  6.4  gives  the  results  after  a  period  of  1  day  and  5  days.  These  data  are 
plotted  in  Eigures  6.5  and  6.6.  As  before,  the  analytic  and  computed 
results  match  very  well. 
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Table  6.4.  Total  head  at  different  elevations  and  after  different  times  for  the  analytic  solution,  SEEP2D, 

SEEP/W,  and  SLIDE  for  the  saturated  flow  problem. 


Total  Head  (ft) 

Time  =  1  day 

Time  =  5  days 

z(ft) 

Analytic 

SEEP2D 

SEEP/W 

SLIDE 

Analytic 

SEEP2D 

SEEP/W 

SLIDE 

48 

72.776 

72.773 

72.775 

72.772 

73.916 

73.915 

73.916 

73.915 

46 

70.580 

70.574 

70.578 

70.572 

72.833 

72.832 

72.832 

72.832 

44 

68.438 

68.430 

68.435 

68.427 

71.752 

71.751 

71.752 

71.751 

42 

66.375 

66.365 

66.371 

66.360 

70.676 

70.674 

70.675 

70.674 

40 

64.411 

64.401 

64.407 

64.395 

69.604 

69.603 

69.604 

69.602 

38 

62.567 

62.556 

62.563 

62.550 

68.539 

68.537 

68.539 

68.537 

36 

60.855 

60.845 

60.851 

60.838 

67.481 

67.479 

67.481 

67.479 

34 

59.287 

59.278 

59.283 

59.270 

66.432 

66.429 

66.431 

66.429 

32 

57.867 

57.857 

57.864 

57.852 

65.391 

65.388 

65.390 

65.388 

30 

56.598 

56.584 

56.596 

56.585 

64.360 

64.357 

64.359 

64.357 

28 

55.478 

55.470 

55.476 

55.467 

63.339 

63.336 

63.339 

63.336 

26 

54.501 

54.498 

54.500 

54.493 

62.328 

62.325 

62.327 

62.325 

24 

53.660 

53.660 

53.660 

53.654 

61.328 

61.325 

61.327 

61.325 

22 

52.944 

52.947 

52.945 

52.941 

60.339 

60.336 

60.338 

60.336 

20 

52.343 

52.351 

52.344 

52.341 

59.360 

59.357 

59.359 

59.357 

Figure  6.5.  Total  head  (ft)  for  different  z  values  (ft)  for  the  analytic,  SEEP/W, 
SLIDE,  and  SEEP2D  solutions  at  Time  =  1  day  for  the  saturated  flow  problem. 
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Figure  6.6.  Total  head  (ft)  for  different  z  values  (ft)  for  the  analytic,  SEEP/W,  SLIDE, 
and  SEEP2D  solutions  at  Time  =  5  days  for  the  saturated  flow  problem. 
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7  Analysis  of  Generic  Levee  Cross  Section 

Before  actual  levee  systems  are  analyzed,  a  generic  levee  cross  section  was 
first  considered  to  provide  a  tool  for  basic  guidance  using  a  simple  exam¬ 
ple  that  contains  the  main  geometrical  elements  of  a  Mississippi  River  and 
Tributaries  (MR&T)  levee  and  also  to  illustrate  the  analysis  process  that 
will  be  used  on  an  actual  levee  analysis. 

7.1  Description  of  the  probiem 

The  geometry,  soil  layers,  hydraulic  conductivities  for  each  layer,  initial 
conditions,  boundary  conditions,  and  the  hydrograph  of  a  river  resulting 
in  a  flood  stage  will  now  be  described. 

7.1.1  Geometry 

Figure  7.1  shows  the  geometry  and  dimensions  of  the  generic  cross  section. 
The  coordinates  of  the  main  points  required  to  define  the  geometry  are 
shown  on  the  figure.  The  cross  section  is  similar  to  a  Lower  Mississippi 
River  Valley  levee  in  that  a  pervious  sand  layer  or  aquifer  is  overlain  by  a 
less  pervious  silty  sand  top  stratum. 

7.1.2  Material  properties 

Table  7.1  gives  the  saturated  hydraulic  conductivities  of  the  three 
materials,  and  Table  7.2  gives  the  unsaturated  and  transient  flow 
parameters.  The  van  Genuchten  (1980)  method  was  used  to  determine  the 
water  content  function  and  the  hydraulic  conductivity  function. 

7.1.3  Initial  and  boundary  conditions 

Both  transient  and  steady-state  analyses  were  performed  using  the  generic 
cross  section.  The  steady-state  boundary  conditions  used  in  SEEP/W  are 
shown  in  Eigure  7.2. 

Eor  the  transient  solution,  the  hydrograph  shown  in  Eigure  7.3  is  applied 
on  the  upstream  boundary  for  a  period  of  28  days.  It  is  important  to  note 
that  the  maximum  elevation  of  the  river  of  17.5  ft  is  reached  in  14  days  and 
then  recedes.  This  type  of  hydrograph  is  similar  to  what  would  be  realized 
with  a  flood  protection  levee. 


Figure  7.1.  Generic  ievee  cross  section. 
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Table  7.1.  Saturated  hydraulic  conductivities. 


Soil  Description 

Horizontal  Saturated  Hydraulic 
Conductivity,  (fl/day) 

Vertical  Saturated  Hydraulic 
Conductivity,  (ft/day) 

Lean  clay 

0.0028 

0.0028 

Silty  sand 

0.28 

0.28 

Uniform  sand 

28. 

28. 

Table  7.2.  Material  properties  of  (1)  volumetric  compressibility  (m^),  (2)  residual  moisture  content 
(Br),  (3)  saturated  moisture  content  (BJ,  and  (4)  van  Genuchten  parameters  (a  and  n). 


Soil  Description 

m,  (ftVIb) 

B^  (ftVftS) 

Bs  (ft3/ft3) 

a(psf) 

n 

Lean  clay 

10-5 

0.05 

0.5 

204.7 

1.23 

Silty  sand 

10-5 

0.035 

0.35 

20.1 

2.24 

Uniform  sand 

10-5 

Not  used 

Not  used 

Not  used 

Not  used 

Figure  7.2.  Steady-state  boundary  conditions. 

Specified  total  head  of  -5  ft 


Figure  7.3.  Hydrograph  where  the  river  elevation  reaches  a 
peak  and  then  decreases. 
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It  is  best  to  have  reasonable  and  verifiable  initial  conditions  when  doing  a 
transient  analysis.  For  simplicity  in  analyzing  the  generic  cross  section,  the 
initial  conditions  are  total  head  =  -5  ft  for  all  nodes  in  the  mesh. 

7.2  Steady-state  analysis 

The  basic  procedure  for  doing  a  transient  seepage  analysis  using  the  generic 
cross  section  as  a  representative  levee  is  described  in  the  following  sections. 
Before  doing  a  transient  analysis,  it  is  recommended  that  a  steady-state 
solution  first  be  done  to  establish  a  correct  mesh  and  boundary  conditions 
and  for  comparing  steady-state  results  with  the  transient  run. 

The  various  programs  use  different  procedures  to  input  the  geometry  data. 
In  SEEP/W  the  procedure  that  is  used  is  (1)  define  points,  (2)  define 
regions,  (3)  define  material  properties,  (4)  define  boundary  conditions, 
and  (5)  define  mesh  characteristics  and  apply  them  to  the  regions. 

Eigure  7.4  shows  a  portion  of  a  mesh  generated  in  SEEP/W  that  is 
primarily  structured  in  areas  other  than  the  levee.  Eurther  refinement  is 
often  used  at  the  toe  of  the  levee. 


Figure  7.4.  Portion  of  the  mostiy  structured  finite  eiement  mesh. 


To  check  boundary  conditions,  it  is  highly  recommended  that  a  run  first  be 
made  where  all  the  saturated  hydraulic  conductivities  are  set  to  the  same 
value.  Eigure  7.5  shows  the  result  of  the  homogeneous  run  with  contours  of 
total  head,  selected  flow  paths,  and  the  phreatic  surface.  It  is  also  advisable 
to  choose  the  optimize/verify  option  in  SEEP/W  before  running  the 
simulation.  The  homogeneous  soil  computation  looked  correct,  so  the  real 
material  properties  were  then  entered  and  the  problem  rerun.  Eigure  7.6 
shows  a  plot  of  total  head  contours  and  selected  flow  paths  for  the  generic 
levee  using  the  specified  soil  properties  given  in  Tables  7.1  and  7.2. 
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Figure  7.5.  Total  head  contours  and  flow  paths  for  the  homogeneous  soil  case. 


Figure  7.6.  Total  head  contours  and  flow  paths  for  the  actual  (anisotropic)  case. 


The  computed  results  from  the  steady-state  finite  element  analysis  are  as 
follows: 

•  Total  head  at  the  toe  (132.0,  o.)  =  0.0  ft. 

•  Total  head  beneath  the  confining  layer  at  the  toe  (132.0,  -10.0)  = 

5.80  ft. 

•  Exit  gradient  at  the  toe  =  0.58. 

•  Flow  rate  through  the  levee  using  the  flux  section  (see  Figure  7.7) 
(132.0,  0.0)  to  (632.0,  0.0)  to  (632.0,  -90.0)  =  65.18  ft3/day/ft. 

SEEP/W,  SEEP2D,  and  SLIDE  provided  essentially  the  same  results. 
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Figure  7.7.  Places  where  data  are  gathered. 


7.3  Transient  analysis 

This  section  presents  primarily  the  results  of  the  transient  analyses  as 
done  by  using  SEEP/W.  A  comparison  of  results  for  the  programs 
SEEP/W,  SLIDE,  and  SEEP2D  is  also  provided,  and  it  is  given  in 
Section  7.3.3. 

7.3.1  Solution  at  different  times 

A  transient  run  using  SEEP/W  was  performed  using  a  time  step  size  of 
0.1  days  and  the  data  provided  in  Tables  7.1  and  7.2.  Table  7.3  shows  the 
results  of  total  head  at  the  toe,  total  head  beneath  the  confining  layer  at  the 
toe,  exit  gradient  at  the  toe,  and  quantity  of  seepage  through  the  same  flux 
section  as  used  in  the  steady-state  analysis  (Eigure  7.7)  for  various  times  for 
the  hydrograph  shown  in  Eigure  7.3.  Eigure  7.7  shows  these  places  of 
measurement.  Eigure  7.8  shows  the  phreatic  surface  after  14  days.  It  is 
important  to  note  that  the  phreatic  surface  is  significantly  different  than  the 
steady-state  position.  The  phreatic  surface  reaches  the  ground  elevation  at 
the  toe  between  14  and  16  days.  Before  that  time,  no  flow  reaches  the 
surface,  so  there  is  no  exit  gradient.  When  the  phreatic  surface  reaches  the 
ground  surface  at  the  toe,  the  total  head  at  the  toe  becomes  zero,  and  the 
exit  gradient  at  the  toe  is  computed  by  dividing  the  total  head  beneath  the 
confining  layer  at  the  toe  by  the  confining  layer’s  thickness  (10  ft).  Table  7.3 
also  shows  a  comparison  with  the  steady-state  solution. 

The  solution  using  the  hydrograph  shown  in  Eigure  7.9  was  run  to 
1,000  days,  and  points  in  the  clay  levee  still  were  not  close  to  achieving 
steady  state.  The  implication  of  this  is  that  no  surface  of  seepage  was 
developed  on  the  downstream  slope  of  the  levee  as  is  seen  in  the  steady- 
state  case.  It  is  important  that  the  engineer  carefully  weigh  all  these  factors 
in  deciding  whether  to  design  based  on  the  conservative  steady-state  solu¬ 
tion  or  use  the  less  conservative  transient  results  in  some  way  in  their 
design. 
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Table  7.3.  Total  head  at  toe  of  levee,  total  head  at  the  bottom  of  the  confining  layer  at  the  toe 
of  levee,  exit  gradient  at  toe  of  levee,  and  flow  rate  (ftVday/ft)  through  the  flux  section  for 
various  times  for  the  hydrograph  In  Figure  7.3. 


Total  Head 
(ft)  at  Toe 

Total  Head  (ft)  Beneath 
Confining  Layer  at  Toe 

Exit  Gradient 
at  Toe 

Flow  Rate 
(ftVday/ft) 

Steady-State 

0.00 

5.80 

0.580 

65.18 

Time  (days) 

River  Elevation  Cycles  Up  and  Down 

0 

-5.00 

-5.00 

0.000 

0.00 

2 

-5.00 

-4.91 

0.000 

0.08 

4 

-5.00 

-4.58 

0.000 

0.75 

6 

-5.00 

-2.73 

0.000 

4.08 

8 

-5.00 

-1.13 

0.000 

8.83 

10 

-5.00 

0.72 

0.000 

13.88 

12 

-4.98 

2.29 

0.000 

19.84 

14 

-1.46 

3.30 

0.000 

28.92 

16 

0.0 

4.26 

0.426 

39.64 

18 

0.0 

4.84 

0.484 

48.01 

20 

0.0 

4.56 

0.456 

50.99 

22 

0.0 

3.29 

0.329 

45.61 

24 

0.0 

1.81 

0.181 

36.54 

26 

0.0 

0.46 

0.046 

27.49 

28 

-0.29 

-0.40 

0.000 

21.90 

Figure  7.8.  Phreatic  surface  after  14  days. 
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Figure  7.9.  Hydrograph  where  the  river  reaches  a  peak  and  remains 

constant. 


7.3.2  Time  to  steady  state 

A  second  hydrograph  is  shown  in  Figure  7.9,  where  the  river  stays  at  its 
maximum  level  of  17.5  ft  from  14  days  and  beyond.  Although  river 
elevation  is  plotted  to  only  30  days,  it  remains  at  17.5  ft  indefinitely.  Using 
this  hydrograph,  the  transient  solution  will  eventually  go  to  a  steady-state 
solution.  Computations  of  time  to  different  percentages  of  steady  state  for 
important  quantities  will  now  be  presented.  Percentages  for  total  heads 
are  computed  for  a  given  node  by 


P 


h^{t)-h^{0) 

Kss-Mo) 


xlOO% 


(7.1) 


where: 


P  =  percentage  of  steady  state 
=  total  head  at  time,  t 
/it(0)  =  total  head  at  time,  t  =  0 
htss  =  total  head  at  steady  state 
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For  flow  rate,  the  equation  used  is 


P 


Qss 


xlOO% 


(7.2) 


where: 


(t)  =  flow  rate  (fts/day/ft)  through  the  flux  section  at  time,  t 
qss  =  steady-state  value  for  flow  rate  (ft3/day/ft)  through  the  flux 
section 


In  the  same  way,  the  equation  used  for  percentages  for  exit  gradient  is 


P 


xlOO% 


i 


ss 


(7.3) 


where: 


i(t)  =  exit  gradient  at  time,  t 
igs  =  steady-state  value  for  exit  gradient  at  the  toe 

Because  the  phreatic  surface  in  the  transient  case  (Figure  7.8)  has  not 
achieved  steady  state,  and  thus  some  of  the  water  entering  the  system  is 
used  for  storage  in  the  unsaturated  zone,  the  time  for  99%  of  steady  state 
for  flow  rate  takes  a  relatively  long  time  to  achieve.  ^  Table  7.4  shows  the 
times  to  a  specific  percentage  of  steady  state  for  the  variables  given  in 
Table  7.3. 


Table  7.4.  Time  to  %  of  steady  state  for  different  percentages  for  quantities. 


Time  to  %  of  Steady  State  (days) 

10% 

20% 

80% 

90% 

95% 

99% 

Total  head  (ft)  at  toe 

13.0 

13.2 

15.0 

15.5 

15.8 

16.0 

Total  head  beneath 
confining  layer  at  toe 

4.7 

5.9 

14.7 

16.9 

19.8 

32.2 

Exit  gradient  at  toe 

15.2 

15.3 

16.7 

19.5 

23.4 

41.0 

Flow  rate  (ftVday/ft) 

6.4 

9.7 

19.1 

22.7 

26.9 

42.5 

1  The  majority  of  the  flow  occurs  through  the  silty  sand  and  sand  layers.  Very  little  flow  occurs  through 
the  levee  embankment. 
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7.3.3  Comparison  of  results  among  the  different  programs 

A  comparison  of  results  between  SEEP/W,  SLIDE,  and  SEEP2D  was  done 
to  see  how  well  they  matched.  To  facilitate  this,  key  points  in  the  levee  were 
selected  as  shown  in  Eigure  7.10.  There  was  a  surprisingly  considerable 
effort  required  to  harmonize  the  different  computer  programs.  Tables  7.5 
through  7.10  with  their  corresponding  plots  in  Eigures  7.11  to  7.16  show  the 
results. 


Figure  7.10.  Selected  points  for  testing  the  different  programs. 


Table  7.5.  Total  head  (ft)  at  different  times  for  the  different  programs  for  point  1. 


Time  (days) 

Total  Head  (ft)  @  Point  1  (73.11,  7.48) 

SEEP2D 

SEEP/W 

SLIDE 

0.0 

-5.0000 

-5.0000 

-5.0000 

0.1 

-5.0000 

-5.0000 

-5.0000 

0.25 

-5.0000 

-5.0000 

-5.0000 

0.5 

-5.0000 

-5.0000 

-5.0000 

1. 

-5.0000 

-5.0000 

-5.0000 

2. 

-5.0000 

-5.0000 

-5.0000 

3. 

-5.0000 

-5.0000 

-5.0000 

4. 

-5.0000 

-5.0000 

-5.0000 

8. 

-5.0000 

-5.0000 

-4.9998 

12. 

-5.0000 

-5.0000 

-4.9759 

16. 

-5.0000 

-5.0000 

-4.9575 

20. 

-5.0000 

-5.0000 

-5.0514 

24. 

-4.9997 

-5.0000 

-5.0803 

28. 

-4.9994 

-5.0000 

-5.0910 
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Figure  7.11.  Total  head  (ft)  at  different  times  for  the  different  programs  for 

point  1. 


Table  7.6.  Total  head  (ft)  at  different  times  for  the  different  programs  for  point  2. 


Time  (days) 

Total  Head  (ft)  @  Point  2  (98.40, 3.53) 

SEEP2D 

SEEP/W 

SLIDE 

0.0 

-5.0000 

-5.0000 

-5.0000 

0.1 

-5.0000 

-5.0000 

-5.0000 

0.25 

-5.0000 

-5.0000 

-5.0000 

0.5 

-5.0000 

-5.0000 

-5.0000 

1. 

-5.0000 

-5.0000 

-5.0000 

2. 

-5.0000 

-5.0000 

-5.0000 

3. 

-5.0000 

-5.0000 

-5.0000 

4. 

-5.0000 

-5.0000 

-5.0000 

8. 

-5.0000 

-5.0000 

-5.0000 

12. 

-4.9998 

-4.9991 

-4.9629 

16. 

-4.9938 

-4.9226 

-4.7663 

20. 

-4.4745 

-4.4533 

-5.0188 

24. 

-4.1293 

-4.0926 

-5.475 

28. 

-4.1110 

-4.0076 

-5.7378 
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Figure  7.12.  Total  head  (ft)  at  different  times  for  the  different  programs  for 

point  2. 
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Table  7.7.  Total  head  (ft)  at  different  times  for  the  different  programs  for  point  3. 


Time  (days) 

Total  Head  (ft)  @  Point  3  (105.60,  -6.67) 

SEEP2D 

SEEP/W 

SLIDE 

0.0 

-5.0000 

-5.0000 

-5.0000 

0.1 

-5.0000 

-5.0000 

-5.0000 

0.25 

-5.0000 

-4.9999 

-5.0000 

0.5 

-5.0000 

-4.9980 

-4.9996 

1. 

-4.9988 

-4.9833 

-4.9939 

2. 

-4.9769 

-4.8944 

-4.9472 

3. 

-4.8511 

-4.7310 

-4.873 

4. 

-4.5703 

-4.4345 

-4.7441 

8. 

-0.9054 

-1.7529 

-1.8772 

12. 

1.4998 

1.7295 

1.2451 

16. 

4.1648 

3.8353 

5.0392 

20. 

5.1057 

5.0894 

5.0388 

24. 

2.2903 

2.1640 

2.0896 

28. 

-0.1384 

-0.1362 

-0.18238 
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Figure  7.13.  Total  head  (ft)  at  different  times  for  the  different  programs  for 

point  3. 


Table  7.8.  Total  head  (ft)  at  different  times  for  the  different  programs  for  point  4. 


Time  (days) 

Total  Head  (ft)  @  Point  4  (132,  -10) 

SEEP2D 

SEEP/W 

SLIDE 

0.0 

-5.0000 

-5.0000 

-5.0000 

0.1 

-5.0000 

-5.0000 

-5.0000 

0.25 

-5.0000 

-4.9998 

-4.9999 

0.5 

-4.9997 

-4.9968 

-4.9992 

1. 

-4.9902 

-4.9756 

-4.9884 

2. 

-4.9099 

-4.8558 

-4.9042 

3. 

-4.4658 

-4.6423 

-4.7841 

4. 

-3.7707 

-4.2367 

-4.5828 

8. 

-0.4260 

-0.8422 

-0.85036 

12. 

2.4355 

2.3660 

2.3684 

16. 

4.1004 

4.1554 

4.4727 

20. 

4.3936 

4.4207 

4.4168 

24. 

1.7240 

1.7111 

1.7192 

28. 

-0.4740 

-0.4552 

-0.42305 
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Figure  7.14.  Total  head  (ft)  at  different  times  for  the  different  programs  for 

point  4. 


Table  7.9.  Total  head  (ft)  at  different  times  for  the  different  programs  for  point  5. 


Time  (days) 

Total  Head  (ft)  @  Point  5  (90.96,  -  21.49) 

SEEP2D 

SEEP/W 

SLIDE 

0.0 

-5.0000 

-5.0000 

-5.0000 

0.1 

-5.0000 

-5.0000 

-5.0000 

0.25 

-5.0000 

-4.9996 

-4.9998 

0.5 

-4.9993 

-4.9951 

-4.995 

1. 

-4.9843 

-4.9654 

-4.9819 

2. 

-4.8765 

-4.8119 

-4.87 

3. 

-4.2801 

-4.5516 

-4.7223 

4. 

-3.4311 

-4.0411 

-4.4626 

8. 

0.4121 

0.0048 

-0.02753 

12. 

3.6684 

3.6025 

3.5706 

16. 

5.4223 

5.4450 

5.7617 

20. 

5.5039 

5.5029 

5.4915 

24. 

2.2752 

2.2475 

2.2576 

28. 

-0.2007 

-0.19168 

-0.17831 
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Figure  7.15.  Total  head  (ft)  at  different  times  for  the  different  programs  for 

point  5. 


Point  5  (90.96,  -21.49) 


Table  7.10.  Total  head  (ft)  at  different  times  for  the  different  programs  for  point  6. 


Time  (days) 

Total  Head  (ft)  @  Point  6  (111.41,  -  21.59) 

SEEP2D 

SEEP/W 

SLIDE 

0.0 

-5.0000 

-5.0000 

-5.0000 

0.1 

-5.00000 

-5.0000 

-5.0000 

0.25 

-5.0000 

-4.9997 

-4.9999 

0.5 

-4.9995 

-4.9960 

-4.999 

1. 

-4.9874 

-4.9707 

-4.9853 

2. 

-4.8936 

-4.8342 

-4.8873 

3. 

-4.3746 

-4.5973 

-4.7532 

4. 

-3.6012 

-4.1400 

-4.523 

8. 

0.0007 

-0.4087 

-0.43166 

12. 

3.0638 

2.9989 

2.9812 

16. 

4.7744 

4.8094 

5.1337 

20. 

4.9710 

4.9832 

4.9755 

24. 

2.0121 

1.9912 

1.9997 

28. 

-0.3345 

-0.3210 

-0.2999 
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Figure  7.16.  Total  head  (ft)  at  different  times  for  the  different  programs  for 

point  6. 


7.3.4  Pore  pressures  in  levee  fill 

Pore  pressures  in  the  levee  fill  are  important  to  slope  stability  computations. 
What  is  shown  next  are  pore  pressures  at  different  times  and  at  different 
points  in  the  levee  fill.  The  percent  of  the  steady-state  value  is  also  given. 
Percentage  of  steady  state  for  pore  pressure  is  computed  by 


P 


u(t)-u(0] 

- ^xl00% 


(7.4) 


where: 


P  =  percentage  of  steady  state 
ii(t)  =  pore  pressure  at  time,  t 
ii(0)  =  pore  pressure  at  time,  t  =  0 
^tss  =  pore  pressure  at  steady  state 

Table  7.11  gives  the  results  for  the  five  points  shown  in  Figure  7.17.  Only 
one  of  the  five  points  achieved  steady  state  in  1,000  days. 
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Table  7.11.  Pore  pressure  (psf)  at  different  times  for  the  different  programs  for  the  original 

soil  parameters. 


Point 

1 

2 

3 

4 

5 

X 

25.9 

66.0 

102.3 

54.2 

82.6 

Z 

2.0 

2.0 

3.9 

6.0 

8.2 

Steady-State 

574.5 

380.2 

172.4 

290.7 

-52.5 

Time  (days) 

Pore  Pressure  (psf)/Percentage  of  Steady-State 

0. 

-436.5 

-436.8 

-555.7 

-684.3 

-824.8 

0.0 

0.0 

0.0 

0.0 

0.0 

50. 

404.9 

131.4 

-517.7 

-684.0 

-824.8 

83.2 

69.6 

5.2 

0.0 

0.0 

100. 

565.1 

249.4 

-62.4 

-661.8 

-824.7 

99.1 

84.0 

67.8 

2.3 

0.0 

150. 

571.6 

299.0 

84.1 

-375.4 

-824.4 

99.7 

90.1 

87.9 

31.7 

0.1 

200. 

572.9 

299.0 

74.6 

-83.6 

-822.1 

99.8 

90.1 

86.6 

61.6 

0.4 

250. 

573.5 

311.5 

82.6 

5.9 

-814.9 

99.9 

91.6 

87.7 

70.8 

1.3 

300. 

573.7 

336.7 

100.4 

66.5 

-796.5 

99.9 

94.7 

0.1 

77.0 

3.7 

350. 

574.0 

343.4 

118.1 

72.2 

-761.0 

99.9 

95.5 

92.5 

77.6 

8.3 

400. 

574.0 

343.6 

133.5 

61.6 

-707.8 

00.0 

95.5 

94.7 

76.5 

15.2 

450. 

574.2 

343.0 

138.9 

91.0 

-646.8 

100.0 

95.4 

95.4 

79.5 

23.0 

500. 

574.2 

342.1 

142.0 

106.9 

-586.4 

100.0 

95.3 

95.8 

81.2 

30.9 

550. 

574.2 

341.1 

145.8 

120.7 

-528.4 

100.0 

95.2 

96.3 

82.6 

38.4 

600. 

574.3 

340.8 

146.9 

144.0 

-472.2 

100.0 

95.2 

96.5 

85.0 

45.7 

650. 

574.3 

341.8 

149.0 

167.3 

-417.5 

100.0 

95.3 

96.8 

87.3 

52.7 

700. 

574.3 

345.2 

150.6 

183.7 

-365.2 

100.0 

95.7 

97.0 

89.0 

59.5 

750. 

574.3 

350.4 

151.6 

187.7 

-365.2 

100.0 

96.4 

97.1 

89.4 

59.5 

800. 

574.3 

354.4 

152.4 

194.7 

-261.5 

100.0 

96.8 

97.2 

90.2 

72.9 

850. 

574.4 

354.4 

154.2 

208.3 

-214.1 

100.0 

96.8 

97.5 

91.6 

79.1 

900. 

574.4 

358.8 

155.7 

222.9 

-173.7 

100.0 

97.4 

97.7 

93.0 

84.3 

950. 

574.4 

360.0 

156.6 

228.7 

-136.3 

100.0 

97.5 

97.8 

93.6 

89.1 

1000. 

574.4 

360.6 

157.3 

229.0 

-111.7 

100.0 

97.6 

97.9 

93.7 

92.3 
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Figure  7.17.  Selected  points  in  the  levee  fill  for  collecting  pore  pressures. 


7.4  Sensitivity  anaiysis 

A  sensitivity  analysis  was  conducted  for  the  critical  values  of  exit  gradient 
at  the  toe,  flow  rate  (see  flux  section  in  Figure  7.7),  and  time  to  steady  state 
using  the  hydrograph  given  in  Figure  7.9. 

7.4.1  Shape  of  the  water  content  function  and  hydraulic  conductivity 
function 

The  shape  and  position  of  the  water  content  function  and  hydraulic 
conductivity  function  depend  on  whether  the  soil  is  becoming  wetter  or 
dryer.  As  shown  in  Figure  4.6,  there  is  a  different  water  content  function  if 
the  soil  is  drying  (desorption)  or  wetting  (adsorption).  The  difference 
between  these  curves  can  be  as  great  as  an  order  of  magnitude  in  matric 
suction  for  the  same  value  of  water  content  (Sleep  2011). 

This  hysteresis  type  behavior  can  be  approximated  by  choosing  a  curve  that 
is  in  the  middle  of  the  two  extremes.  For  this  study,  it  was  assumed  that  the 
van  Genuchten  (1980)  procedure  provided  a  curve  that  could  be  considered 
as  an  average  between  the  drying  and  wetting  curves.  Figures  7.18  and  7.19 
show  the  curves  used  for  this  analysis  where  only  the  levee  fill  data  were 
modified.  It  is  important  to  note  that  only  the  clay  material  where  the 
position  of  the  phreatic  surface  changes  significantly  is  considered. 

To  determine  the  curves,  the  van  Genuchten  (1980)  formulation  is  used.  It 
is  given  by 


a -a 


1  + 


0  =  0,4 


m 


(7.5) 
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Figure  7.18.  Volumetric  moisture  content  for  the  levee  fill  as  a  function  of 
suction  pressure  for  the  wetting,  original,  and  drying  cases. 


Figure  7.19.  Hydraulic  conductivity  function  for  the  levee  fill  as  a  function  of 
suction  pressure  for  the  wetting,  original,  and  drying  cases. 
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where  from  Table  7.2,  9^  =  0.05,  Og  =  0.5,  a  =  204.7  psf,  and  n  =1.23.  Also, 

m  =  l--  (7.6) 

n 

The  xjj  value  (xpave)  at  the  average  moisture  content  value  (Gave)  between  9^. 
and  9s  was  determined.  That  is, 

(7-7) 

Equations  7.5  and  7.7  are  used  to  compute  xpave  the  following  alge¬ 
braic  operations: 


^ave  —  + 


0 


1  + 


Wa 


a 


|(0,+e,)-e,= 


1  + 


Wa 


a 


1  + 


ave 

a 


1+ 

y^ave 

a 

1+ 

a 

=  2 


=  2 


w  =a 

T  ave 


2  m  _1 


(7.8) 


For  this  study,  the  difference  in  xp  between  the  wetting  and  drying  curves 
at  9 ave  is  assumed  to  be  approximately  an  order  of  magnitude  (Figure  7.18 
at  0  =  0.275).  Values  of  ip  one-half  order  of  magnitude  on  either  side  of 
xpave  were  next  computed  using 
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1 

logioV^±=logioV^a.e±2 

Making  both  sides  of  Equation  7.9  a  power  of  10  yields 

1  1 
j_QlogioV'+  ^  j_QlogioV'-  ^ 

This  yields  the  result, 

ip^  =  (7.10) 

w_  =  ^^  (7.11) 

Vio 

Finally,  new  values  of  a  in  Equation  7.5  are  computed  so  that  9  -  Gave  at 
xp  =  xIj+.  This  is  accomplished  by 

Vave  _  y^+ 

a  a 


This  gives 


a^=a-^  (7.12) 

Wave 

The  hydraulic  conductivity  function  was  computed  from  the  van 
Genuchten  (1980)  function.  Starting  with  the  definition  of  effective 
saturation,  Sg, 


S  = 


O-Or 

0 


the  hydraulic  conductivity  function  for  x(j  >0  is  given  by 


K=K4^e 


1- 


1-s: 


(7.13) 


(7.14) 
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where 


kg  =  saturated  hydraulic  conductivity 

Values  for  this  particular  problem  are  xpave  —  4084.8  psf,  a_  =  647.3  psf, 
a  =  204.7  psf,  and  a+  =  64.7  psf .  Figures  7.14  and  7.15  give  the  resulting 
curves  for  the  water  content  and  hydraulic  conductivity  functions  for  the 
above  values  of  a. 

Table  7.12  gives  the  exit  gradient  at  the  toe  of  the  levee  at  different  times 
and  for  the  wetting,  original,  and  drying  curves  for  the  hydrograph  that 
goes  up  and  remains,  and  Table  7.13  gives  the  same  type  data  for  flow  rate 
through  the  flux  section  given  in  Figure  7.7.  Tables  7.14  and  7.15  give  the 
time  to  steady  state  for  the  wetting,  original,  and  drying  curves. 

Figures  7.20  through  7.23  show  plots  of  the  results  found  in  the  tables.  For 
these  key  design  parameters,  the  impact  of  differences  in  unsaturated  flow 
curves  seems  small.  This  is  reasonable  because  the  levee  fill  material  has  a 
very  low  permeability,  and  little  flow  occurs  through  this  material. 

The  results  for  the  points  in  Figure  7.17  were  tabulated  as  in  Table  7.11 
from  the  wetting  and  drying  computations  discussed  above.  Tables  7.16 
and  7.17  give  the  results  of  this  tabulation,  and  Figures  7.24  to  7.28  show 
plots  of  wetting,  original,  and  drying  curves  for  the  five  points  in 
Figure  7.17.  From  these  results,  it  is  clear  that  the  hydraulic  conductivity 
and  volumetric  water  content  functions  for  the  levee  fill  make  a  significant 
impact  on  the  results  of  pore  pressures  in  the  drier  part  of  the  levee  fill. 


Table  7.12.  Exit  gradient  at  the  toe  of  the  levee  at  different  times  and  for  wetting,  original, 
and  drying  curves  for  the  hydrograph  that  goes  up  and  remains. 


Time  (days) 

Exit  Gradient  at  Toe 

Wetting 

Original 

Drying 

0 

0.000 

0.000 

0.000 

2 

0.000 

0.000 

0.000 

4 

0.000 

0.000 

0.000 

6 

0.000 

0.000 

0.000 

8 

0.000 

0.000 

0.000 

10 

0.000 

0.000 

0.000 

12 

0.000 

0.000 

0.000 

14 

0.000 

0.000 

0.000 
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Time  (days) 

Exit  Gradient  at  Toe 

Wetting 

Original 

Drying 

16 

0.446 

0.433 

0.416 

18 

0.504 

0.500 

0.494 

20 

0.530 

0.528 

0.526 

22 

0.545 

0.544 

0.543 

24 

0.555 

0.554 

0.553 

26 

0.561 

0.560 

0.559 

28 

0.566 

0.564 

0.564 

30 

0.569 

0.567 

0567 

35 

0.573 

0.571 

0.571 

40 

0.576 

0.574 

0.573 

45 

0.577 

0.575 

0.575 

50 

0.578 

0.576 

0.575 

55 

0.578 

0.577 

0.575 

60 

0.579 

0.577 

0.576 

Table  7.13.  Flow  rate  (ftVday/ft)  through  the  flux  section  at  different  times  and  for  wetting, 
original,  and  drying  curves  for  the  hydrograph  that  goes  up  and  remains. 


Time  (days) 

Flow  Rate  (ftyday/ft) 

Wetting 

Original 

Drying 

0 

0.00 

0.00 

0.00 

2 

0.08 

0.08 

0.08 

4 

0.75 

0.75 

0.75 

6 

4.08 

4.08 

4.08 

8 

8.83 

8.83 

8.83 

10 

13.89 

13.89 

13.88 

12 

19.86 

19.84 

19.82 

14 

29.15 

28.93 

28.83 

16 

40.57 

39.81 

39.06 

18 

49.19 

48.78 

48.18 

20 

54.53 

54.34 

54.00 

22 

57.92 

57.74 

57.58 

24 

60.08 

60.01 

59.89 

26 

61.55 

61.45 

61.39 
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Time  (days) 

Flow  Rate  (ftVday/ft) 

Wetting 

Original 

Drying 

28 

62.52 

62.38 

62.37 

30 

63.19 

63.06 

63.03 

35 

64.14 

63.96 

63.97 

40 

64.56 

64.39 

64.39 

45 

64.76 

64.64 

64.58 

50 

64.88 

64.77 

64.64 

55 

64.95 

64.85 

64.71 

60 

65.00 

64.90 

64.78 

Table  7.14.  Time  to  steady  state  (days)  for  exit  gradient  at  the  toe  of  the  levee  for  wetting, 

original,  and  drying  curves. 


%  of  Steady  State 

Time  to  %  of  Steady  State  for  Exit  Gradient  at  the  Toe  (days) 

Wetting 

Original 

Drying 

10 

15.1 

15.2 

15.1 

20 

15.3 

15.3 

15.2 

80 

16.5 

16.7 

17.1 

90 

19.3 

19.5 

19.7 

95 

23.2 

23.4 

23.5 

99 

36.4 

41.0 

42.0 

Table  7.15.  Time  to  steady  state  (days)  for  flow  rate  (ftVday/ft)  through  the  flux  section  for 

wetting,  original,  and  drying  curves. 


%  of  Steady  State 

Time  to  %  of  Steady  State  for  Flow  Rate  (days) 

Wetting 

Original 

Drying 

10 

6.5 

6.4 

7.0 

20 

9.7 

9.7 

9.7 

80 

19.0 

19.1 

19.3 

90 

22.6 

22.7 

22.8 

95 

26.7 

26.9 

27.0 

99 

38.5 

42.5 

43.5 

ERDC  TR-16-8 


59 


Figure  7.20.  Exit  gradient  at  the  toe  of  the  ievee  at  different  times  and  for  wetting, 
originai,  and  drying  curves  for  the  hydrograph  that  goes  up  and  remains. 


Figure  7.21.  Fiow  rate  (ftVday/ft)  through  the  fiux  section  at  different  times  and  for 
wetting,  originai,  and  drying  curves  for  the  hydrograph  that  rises  and  remains  constant. 
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Figure  1.22.  Time  to  steady  state  (days)  for  exit  gradient  at  the  toe  of  the  ievee  for 
wetting,  originai,  and  drying  curves. 
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Figure  7.23.  Time  to  steady  state  (days)  for  fiow  rate  (ftVday/ft)  through  the  fiux 
section  for  wetting,  originai,  and  drying  curves. 
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Table  7.16.  Pore  pressure  (psf)  at  different  times  for  the  different  programs  for  the  wetting 

soil  property  curves  for  the  levee  fill. 


Point 

1 

2 

3 

4 

5 

X 

25.9 

66.0 

102.3 

54.2 

82.6 

z 

2.0 

2.0 

3.9 

6.0 

8.2 

Steady-State 

574.5 

382.8 

172.5 

295.3 

-50.7 

Time  (days) 

Pore  Pressure  (psf)/Percentage  of  Steady-State 

0. 

-436.5 

-436.8 

-555.7 

-684.3 

-824.8 

0.0 

0.0 

0.0 

0.0 

0.0 

50. 

562.6 

228.9 

-96.5 

-634.8 

-824.5 

98.8 

81.2 

63.1 

5.0 

0.0 

100. 

572.6 

274.9 

60.6 

-104.9 

-805.6 

99.8 

86.8 

84.6 

59.1 

2.5 

150. 

573.6 

324.2 

94.0 

31.6 

-709.9 

99.9 

92.8 

89.2 

73.1 

14.8 

200. 

573.9 

330.4 

118.5 

55.0 

-545.7 

99.9 

93.6 

92.6 

75.5 

36.1 

250. 

574.1 

334.5 

130.2 

97.0 

-400.0 

100.0 

94.1 

94.2 

79.8 

54.9 

300. 

574.2 

338.2 

138.3 

133.0 

-298.1 

100.0 

94.6 

95.3 

83.4 

68.0 

350. 

574.3 

342.2 

144.9 

162.7 

-237.7 

100.0 

95.1 

96.2 

86.5 

75.8 

400. 

574.3 

345.7 

150.0 

189.8 

-203.2 

100.0 

95.5 

96.9 

89.2 

80.3 

450. 

574.4 

347.8 

153.7 

202.9 

-182.0 

100.0 

95.7 

97.4 

90.6 

83.0 

500. 

574.4 

349.6 

156.8 

211.6 

-166.2 

100.0 

96.0 

97.8 

91.5 

85.1 

550. 

574.4 

352.0 

159.1 

226.4 

-154.3 

100.0 

96.2 

98.2 

93.0 

86.6 

600. 

574.4 

354.9 

161.0 

237.8 

-144.8 

100.0 

96.6 

98.4 

94.1 

87.8 

650. 

574.4 

356.7 

162.5 

241.3 

-135.7 

100.0 

96.8 

98.6 

94.5 

89.0 

700. 

574.4 

356.7 

163.6 

242.8 

-125.9 

100.0 

96.8 

98.8 

94.6 

90.3 

750. 

574.4 

358.9 

164.5 

243.8 

-122.2 

100.0 

97.1 

98.9 

94.7 

90.8 

800. 

574.4 

359.7 

165.3 

246.8 

-120.1 

100.0 

97.2 

99.0 

95.1 

91.0 

850. 

574.4 

360.5 

165.9 

250.1 

-117.1 

100.0 

97.3 

99.1 

95.4 

91.4 

900. 

574.5 

361.3 

166.6 

252.6 

-113.1 

100.0 

97.4 

99.2 

95.6 

91.9 

950. 

574.5 

361.9 

167.1 

254.5 

-108.5 

100.0 

97.4 

99.3 

95.8 

92.5 

1000. 

574.5 

362.4 

167.5 

256.8 

-104.3 

100.0 

97.5 

99.3 

96.1 

93.1 
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Table  7.17.  Pore  pressure  (psf)  at  different  times  for  the  different  programs  for  the  drying  soil 

property  curves  for  the  levee  fill. 


Point 

1 

2 

3 

4 

5 

X 

25.9 

66.0 

102.3 

54.2 

82.6 

z 

2.0 

2.0 

3.9 

6.0 

8.2 

Steady-State 

574.5 

377.4 

172.5 

283.9 

-53.0 

Time  (days) 

Pore  Pressure  (psf)/Percentage  of  Steady-State 

0. 

-436.5 

0.0 

-436.8 

0.0 

-555.7 

0.0 

-684.3 

0.0 

-824.8 

0.0 

50. 

167.4 

59.7 

-24.4 

50.6 

-555.3 

0.0 

-684.3 

-0.0 

-824.8 

0.0 

100. 

530.5 

95.6 

205.2 

78.8 

-289.5 

36.6 

-684.0 

0.0 

-824.8 

0.0 

150. 

563.1 

98.9 

236.5 

82.7 

-85.5 

64.6 

-682.1 

0.2 

-824.8 

0.0 

200. 

570.6 

99.6 

302.7 

90.8 

69.4 

85.8 

-655.0 

3.0 

-824.8 

0.0 

250. 

572.1 

99.8 

312.5 

92.0 

110.8 

91.5 

-304.3 

39.2 

-824.8 

0.0 

300. 

572.9 

99.8 

279.8 

88.0 

114.6 

92.0 

-79.9 

62.4 

-824.8 

0.0 

350. 

573.3 

99.9 

299.7 

90.5 

78.4 

87.1 

-23.4 

68.3 

-824.8 

0.0 

400. 

573.6 

99.9 

313.8 

92.2 

93.9 

89.2 

21.3 

72.9 

-824.2 

0.1 

450. 

573.8 

99.9 

330.6 

94.3 

102.1 

90.3 

69.0 

77.8 

-823.3 

0.2 

500. 

574.0 

100.0 

347.2 

96.3 

110.3 

91.5 

105.7 

81.6 

-821.5 

0.4 

550. 

574.0 

100.0 

353.6 

97.1 

123.5 

93.3 

116.2 

82.7 

-818.7 

0.8 

600. 

574.1 

100.0 

356.4 

97.4 

134.8 

94.8 

73.4 

78.3 

-814.0 

1.4 

650. 

574.2 

100.0 

357.6 

97.6 

145.6 

96.3 

75.9 

78.5 

-807.4 

2.3 

700. 

574.2 

100.0 

358.3 

97.7 

149.0 

96.8 

104.5 

81.5 

-799.0 

3.3 

750. 

574.3 

100.0 

358.7 

97.7 

153.2 

97.3 

116.2 

82.7 

-789.6 

4.6 

800. 

574.2 

100.0 

358.6 

97.7 

155.8 

97.7 

118.2 

82.9 

-779.7 

5.8 

850. 

574.3 

100.0 

358.4 

97.7 

155.8 

97.7 

133.4 

84.5 

-769.7 

7.1 

900. 

574.3 

100.0 

357.9 

97.6 

158.6 

98.1 

149.2 

86.1 

-759.5 

8.5 

950. 

574.3 

100.0 

357.3 

97.5 

159.9 

98.3 

153.5 

86.5 

-749.2 

9.8 

1000. 

574.3 

100.0 

356.7 

97.5 

159.6 

98.2 

163.3 

87.5 

-739.0 

11.1 
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Figure  7.24.  Pore  pressure  (psf)  at  Point  1  at  different  times. 


Figure  7.25.  Pore  pressure  (psf)  at  Point  2  at  different  times. 
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Figure  7.26.  Pore  pressure  (psf)  at  Point  3  at  different  times. 


Figure  7.27.  Pore  pressure  (psf)  at  Point  4  at  different  times. 
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Figure  7.28.  Pore  pressure  (psf)  at  Point  5  at  different  times. 


7.4.2  Value  of  volumetric  compressibility 

The  volumetric  compressibility,  m^,  plays  an  important  role  in  the  solution 
of  transient  seepage  in  SEEP/W  as  it  is  on  the  right-hand  side  of  the 
governing  equation  for  transient  flow  (Equation  3.5).  It  is  used  in  SEEP/W 
in  a  region  that  has  saturated  flow.  The  value  of  was  varied  by  four 
orders  of  magnitude  by  using  values  of  10“^,  10“^,  and  10“^  ft^/lb.  The 
values  of  volumetric  compressibility  correspond  to  very  soft  clays  or  peats, 
plastic  clays,  and  dense  sands,  respectively.  The  exit  gradient  at  the  toe 
and  flow  rate  through  the  flux  section  given  in  Eigure  7.7  for  different 
times  were  calculated  and  plotted.  Tables  7.18  through  7.21  and 
Eigures  7.29  through  7.32  show  the  results. 

A  significant  difference  was  observed  for  compressibility  values  ranging 
from  10“^  ft^/lb  to  10“^  ft^/lb,  but  very  little  in  comparison  from 
10“^  ft^/lb  to  10“^  ft^/lb.  This  can  be  understood  by  examining  the 
solution  to  the  i-D  problem  shown  in  Eigure  6.1  where  a  column  of  soil  of 
length,  L,  saturated  hydraulic  conductivity,  kg,  density,  Yw^  and  volumetric 
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compressibility,  ruy,  has  an  initial  total  head  of  Hq  until  a  head  of  Hi  is 
applied  at  the  top  such  that  the  increased  head  values  propagate  down¬ 
ward  as  time  increases.  The  solution  as  given  in  Equation  B.22  is 


A,=«„+ 


k=l 


z  +  5];(-lf  —  sin(A,z)e 

A,, 


(7.14) 


where 


=  =  0,1,2, 


Yu,m 


Table  7.18.  Exit  gradient  at  the  toe  of  the  levee  at  different  times  and  for  different  volumetric 
compressibility  values  for  the  hydrograph  that  goes  up  and  remains. 


Time  (days) 

Exit  Gradient  at  Toe 

rUy  =  10“^  /lb 

rriy  =  10  ^ 

rriy  =  10  ^  /lb 

0 

0.000 

0.000 

0.000 

2 

0.000 

0.000 

0.000 

4 

0.000 

0.000 

0.000 

6 

0.000 

0.000 

0.000 

8 

0.000 

0.000 

0.000 

10 

0.000 

0.000 

0.000 

12 

0.000 

0.000 

0.000 

14 

0.000 

0.000 

0.510 

16 

0.000 

0.433 

0.537 

18 

0.000 

0.500 

0.550 

20 

0.000 

0.528 

0.558 

22 

0.000 

0.544 

0.563 

24 

0.000 

0.554 

0.566 

26 

0.000 

0.560 

0.569 

28 

0.000 

0.564 

0.571 

30 

0.000 

0.567 

0.571 

35 

0.000 

0.571 

0574 

40 

0.000 

0.574 

0.575 

45 

0.000 

0.575 

0.576 

50 

0.000 

0.576 

0.577 

55 

0.000 

0.577 

0.577 

60 

0.000 

0.577 

0.578 

80 

0.000 

0.578 

0.578 
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Time  (days) 

Exit  Gradient  at  Toe 

=  10“^  /lb 

=  10“^ 

=  10  ^  /lb 

100 

0.000 

0.579 

0.579 

150 

0.000 

0.579 

0.580 

200 

0.095 

0.580 

0.580 

300 

0.321 

0.580 

0.580 

400 

0.444 

0.580 

0.580 

600 

0.544 

0.580 

0.580 

800 

0.570 

0.580 

0.580 

1000 

0.577 

0.580 

0.580 

Figure  7.29.  Exit  gradient  at  the  toe  of  the  ievee  at  different  times  and  for  different  voiumetric 
compressibiiity  vaiues  for  the  hydrograph  that  rises  and  remains  constant. 
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Table  7.19.  Flow  rate  (ftVday/ft)  through  the  flux  section  at  different  times  and  for  different 
volumetric  compressibility  values  for  the  hydrograph  that  rises  and  remains  constant. 


Time  (days) 

Flow  Rate  (ftyday/ft) 

=  10  ^  jib 

=  10  ^  ft^ fib 

=  10  ’  fib 

0 

0.00 

0.00 

0.00 

2 

0.00 

0.08 

0.99 

4 

0.00 

0.75 

3.80 

6 

0.00 

4.08 

10.17 

8 

0.00 

8.83 

15.82 

10 

0.00 

13.89 

23.04 

12 

0.00 

19.84 

34.29 

14 

0.00 

28.93 

47.85 

16 

0.00 

39.81 

54.74 

18 

0.00 

48.78 

58.45 

20 

0.00 

54.34 

60.66 

22 

0.00 

57.74 

61.99 

24 

0.00 

60.01 

62.86 

26 

0.00 

61.45 

63.43 

28 

0.00 

62.38 

63.81 

30 

0.00 

63.06 

64.04 

35 

0.00 

63.96 

64.44 

40 

0.00 

64.39 

64.67 

45 

0.01 

64.64 

64.79 

50 

0.02 

64.77 

64.87 

55 

0.05 

64.85 

64.91 

60 

0.08 

64.90 

64.94 

80 

0.44 

65.00 

65.02 

100 

1.24 

65.04 

65.06 

150 

5.67 

65.11 

65.11 

200 

12.77 

65.11 

65.11 

300 

33.46 

65.14 

65.14 

400 

46.77 

65.15 

65.15 

600 

60.27 

65.16 

65.16 

800 

63.86 

65.16 

65.17 

1000 

64.79 

65.17 

65.17 
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Figure  7.30.  Flow  rate  (ft3/day/ft)  through  the  flux  section  at  different  times  and  for  different 
volumetric  compressibility  values  for  the  hydrograph  that  rises  and  remains  constant. 


Table  7.20.  Time  to  steady  state  (days)  for  exit  gradient  at  the  toe  of  the  levee  for  different 

values  of  volumetric  compressibility. 


%  of  steady 

State  Exit 
Gradient 

Time  (days) 

=  10  ^  jib 

=  10  ^  /lb 

=  10  ^  ft^/ lb 

10 

188.3 

15.2 

13.1 

20 

205.4 

15.3 

13.2 

80 

419.9 

16.7 

13.9 

90 

512.5 

19.5 

14.3 

95 

605.0 

23.4 

17.0 

99 

762.0 

41.0 

25.0 
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Figure  7.31.  Time  to  steady  state  (days)  for  exit  gradient  at  the  toe  of  the  ievee  for  different 

vaiues  of  voiumetric  compressibiiity. 


Tabie  7.21.  Time  to  steady  state  (days)  for  fiow  rate  through  the  fiux  section  for  different 

vaiues  of  voiumetric  compressibiiity. 


%  of  steady 
State  Fiow  Rate 

Time  (days) 

=  10  ^  ft^  jib 

=  10  ^  /lb 

=  10  ’  ft^ /lb 

10 

156.7 

6.4 

4.7 

20 

201.5 

9.7 

7.1 

80 

453.2 

19.1 

15.1 

90 

557.7 

22.7 

18.2 

95 

661.5 

26.9 

21.9 

99 

914.5 

42.5 

36.0 
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Figure  7.32.  Time  to  steady  state  (days)  for  fiow  rate  (ft3/day/ft)  through  the  fiux  section  for 
different  vaiues  of  voiumetric  compressibiiity. 


Since  Equation  7.14  typically  converges  very  quickly,  it  is  acceptable  to 
examine  the  result  of  only  using  the  first  term  of  the  series.  Thus, 


H  1  ^^-^0 

'  tin  H - 1 - 


jT^k 

2V 

JIZ 

— 2 — 

z  — 

sin 

h 

Now  suppose  that  everything  is  held  the  same  except  ruy.  Then  to  get  the 
same  value  of  for  say  rriy^,  niy^,  and  m^y,  at  times  t3,  ts,  and  ty, 
respectively. 


m 


l?3 


m 


v5 


m 


v7 
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Thus,  if  171^2  =  10“^  ftVlb,  =  10“^  ftVlb,  and  m^-j  =  10“’^  ftVlb,  the 
times  to  reach  that  same  total  head  are  =  lOOOOty  and  —  lOOty.  So  for 
this  i-D  problem,  plotting  times  to  percentage  of  steady  state  on  a  linear 
scale  shows  an  illusion  that  the  result  for  10“^  ft^/lb  is  much  different 
from  the  results  for  10“^  and  10“^  ft^/lb.  However,  since  is  varied  on  a 
log  scale,  the  results  would  appear  more  similar  if  a  log  scale  was  also  used 
for  the  time  axis.  Although  the  generic  levee  is  more  complicated,  the  same 
phenomenon  is  exhibited. 

7.4.3  Value  of  saturated  hydraulic  conductivity  of  the  confining  layer 

Saturated  hydraulic  conductivity  for  the  confining  layer  was  varied  so  that 
50%  less  than,  the  same  as,  and  50%  greater  than  its  original  value  were 
used  in  the  computer  runs.  The  exit  gradients  at  the  toe  and  flow  rates  on 
the  land  side  at  different  times  was  tabulated  and  plotted  as  before. 

Tables  7.22  and  7.23  and  Figures  7.33  and  7.34  show  the  results. 


Table  7.22.  Exit  gradient  at  the  toe  of  the  levee  at  different  times  and  for  different  hydraulic 
conductivity  values  for  the  hydrograph  that  goes  up  and  remains. 


Exit  Gradient 

Time  (days) 

fej.  =  0. 14  ft/ day 

fej.  =  0. 28  ft/day 

fej.  =  0. 42  ft/day 

0 

0.000 

0.000 

0.000 

2 

0.000 

0.000 

0.000 

4 

0.000 

0.000 

0.000 

6 

0.000 

0.000 

0.000 

8 

0.000 

0.000 

0.000 

10 

0.000 

0.000 

0.000 

12 

0.000 

0.000 

0.000 

14 

0.000 

0.000 

0.430 

16 

0.000 

0.433 

0.504 

18 

0.000 

0.500 

0.534 

20 

0.000 

0.528 

0.549 

22 

0.000 

0.544 

0.558 

24 

0.492 

0.554 

0.564 

26 

0.514 

0.560 

0.567 

28 

0.527 

0.564 

0.569 

30 

0.534 

0.567 

0.571 

35 

0.546 

0.571 

0.573 
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Exit  Gradient 

Time  (days) 

=  0. 14  ft/day 

=  0. 28  ft/day 

k^  =  0. 42  ft/day 

40 

0.551 

0.574 

0.575 

45 

0.554 

0.575 

0.576 

50 

0.557 

0.576 

0.576 

55 

0.559 

0.577 

0.576 

60 

0.560 

0.577 

0.577 

Figure  7.33.  Exit  gradient  at  the  toe  of  the  ievee  at  different  times  and  for  different  hydrauiic 
conductivity  vaiues  of  the  confining  iayer  for  the  hydrograph  that  goes  up  and  remains. 
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Table  7.23.  Flow  rate  (ftVday/ft)  through  the  flux  section  at  different  times  and  for  different 
hydraulic  conductivity  values  for  the  hydrograph  that  goes  up  and  remains. 


Time  (days) 

Flow  Rate  (ftyday/ft) 

=  0. 14  ft/day 

=  0. 28  ft/day 

fej.  =  0.42  ft/day 

0 

0.00 

0.00 

0.00 

2 

0.12 

0.08 

0.05 

4 

1.33 

0.75 

0.49 

6 

3.88 

4.08 

3.98 

8 

9.87 

8.83 

8.06 

10 

15.53 

13.89 

13.15 

12 

21.47 

19.84 

24.05 

14 

26.04 

28.93 

40.66 

16 

29.32 

39.81 

53.33 

18 

33.00 

48.78 

60.61 

20 

37.94 

54.34 

64.86 

22 

41.81 

57.74 

67.40 

24 

45.46 

60.01 

68.99 

26 

48.37 

61.45 

69.97 

28 

50.35 

62.38 

70.62 

30 

51.68 

63.06 

71.00 

35 

53.65 

63.96 

71.54 

40 

54.62 

64.39 

71.83 

45 

55.12 

64.64 

71.97 

50 

55.45 

64.77 

72.05 

55 

55.66 

64.85 

72.10 

60 

55.80 

64.90 

72.14 
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Figure  7.34.  Flow  rate  (ftVday/ft)  exiting  the  levee  on  the  land  side  at  different  times  and  for 
different  hydraulic  conductivity  values  of  the  confining  layer  for  the  hydrograph  that  goes  up 

and  remains. 
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8  Use  of  Transient  Seepage  Analyses  in 
Practice 

There  are  several  scenarios  when  a  transient  seepage  analysis  may  be 
worth  the  additional  effort  to  conduct.  The  pore  pressures,  flow  rates,  and 
gradients  determined  from  a  transient  analysis  will  be  the  most  accurate 
and  reliable  for  cross  sections  that  have  these  characteristics: 

•  Sandy,  as  opposed  to  clayey,  soils 

•  Homogeneous  soils  with  well-defined  soil  properties 

•  Accurate  pre-event  hydrograph 

•  Well-known  initial  conditions  (initial  water  content  and  matric 
suction) 

•  Well-defined  hydraulic  boundary  conditions  (initial  and  final) 

•  Areas  of  interest  remotely  situated  from  areas  of  stress  induced  pore 
pressure  changes 

•  Previous  calibration  of  pore  pressures  and  gradient  to  hydrographs  and 
soil  properties  from  installed  instrumentation 

The  engineer  must  decide  if  the  amount  of  information  that  exists  at  the 
site  is  adequate  to  provide  an  answer  that  is  suitable  for  their  needs.  As 
more  and  more  transient  analyses  are  conducted  on  USAGE  projects,  the 
knowledge  base  will  increase  and  the  analyses  should  become  more 
reliable. 

8.1  Determination  of  time  to  steady  state 

It  is  difficult  to  determine  a  specific  time  at  which  a  transient  flow 
situation  achieves  a  steady  state.  In  real  world  problems,  the  hydraulic 
boundary  conditions,  defined  by  the  hydrograph,  rarely  become  constant 
over  time. 

As  shown  in  the  previous  example,  the  relative  percentage  of  a  steady-state 
quantity  varies  spatially  within  the  levee.  In  the  example,  the  exit  hydraulic 
gradient,  and  the  associated  factor  of  safety  against  uplift,  reached  about 
95%  of  the  steady-state  value  in  25  days.  However,  it  takes  much  longer  for 
the  position  of  the  phreatic  surface  within  the  levee  to  achieve  a  steady-state 
condition.  The  factor  of  safety  for  stability  of  the  landslide  slope  decreases 
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with  increasing  elevation  of  the  phreatic  surface,  and  it  would  take  more 
than  1,000  days  before  the  factor  of  safety  reaches  the  minimum  value, 
assuming  that  there  are  no  shear-induced  pore  pressures  present. 

The  estimates  of  the  time  for  steady  state  should  be  considered  only 
approximate,  and  these  calculations  should  be  performed  using  conserva¬ 
tive  assumptions  when  possible.  As  an  example,  for  the  calculation  of  exit 
gradient  for  the  generic  levee,  it  would  be  conservative  to  assume  that  the 
phreatic  surface  on  the  land  side  is  at  the  ground  surface.  Making  this 
assumption  would  give  the  engineer  more  confidence  that  the  time  to 
steady  state  and  exit  gradient  would  not  be  underestimated. 

8.2  Use  and  misuse  of  transient  anaiyses  in  siope  stabiiity 

There  has  been  a  growing  trend  to  use  transient  analyses  to  calculate  pore 
pressures  that  are  later  used  in  effective  stress  stability  analyses.  These 
uncoupled  analyses  are  normally  an  incorrect  procedure  and  can  lead  to 
misleading  results.  Unfortunately,  this  method  has  been  recommended  as  a 
viable  analysis  method  in  a  variety  of  published  sources  (ILIT  2006;  Geo- 
Slope  2012;  Fredlund  et  al.  2011;  and  USER  2011).  This  procedure  ignores 
the  pore  pressures  that  are  generated  by  both  changes  in  total  stress  and 
shear  stress  and,  depending  on  the  nature  of  the  soil,  these  methods  may  be 
conservative  or  non-conservative.  Details  regarding  this  error  can  be  found 
in  VandenBerge  et  al.  (2015).  The  different  mechanisms  of  pore  pressure 
development  in  dams  and  levees  upon  hydraulic  loading  can  be  explained 
using  the  example  in  Figures  8.1  and  8.2. 

Figure  8.1a  shows  a  levee  with  an  impermeable  membrane  located  on  the 
riverside  slope.  Figure  8.1b  shows  a  flood  loading  on  the  slope.  If  the  levee 
fill  material  is  relatively  saturated  prior  to  the  flood  loading,  it  would  be 
expected  that  the  pore  pressure  would  increase  at  Point  A  due  to  the 
boundary  water  pressure  applied  to  the  membrane.  The  pore  pressure 
would  increase  rapidly  to  the  maximum  value  and  then  dissipate  over 
time,  as  shown  in  Figure  8.1c.  A  similar  example  is  given  in  Bishop  and 
Bjerrum  (i960).  This  change  in  pore  pressure  is  solely  due  to  the  change  in 
stress  because  there  is  no  flow,  owing  to  the  presence  of  the  impermeable 
membrane. 
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Figure  8.1.  Pore  pressure  development  and  dissipation  due  solely  to  change  in 

boundary  pressures. 


Au 


Time 

(c)  Change  in  pore  pressure  due  to  boundary  water  pressures. 


Figure  8.2a  shows  the  same  cross  section,  but  for  this  case,  there  is  no 
impermeable  membrane.  After  the  water  level  rises,  a  seepage  front  would 
be  initiated,  and  the  pore  pressure  would  gradually  increase  at  Point  A. 
The  pore  pressure  change  over  time,  due  solely  to  the  change  in  hydraulic 
boundary  conditions,  is  shown  in  Figure  8.2b.  After  enough  time  has 
passed,  the  pore  pressure  would  reach  the  maximum  value,  which  would 
be  equivalent  to  the  pore  pressure  calculated  from  a  steady-state  seepage 
analysis.  The  time  required  to  reach  the  maximum  pore  pressure  may  be 
days  to  years,  depending  on  the  properties  of  the  levee  fill  material. 
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Figure  8.2.  Pore  pressure  development  and  dissipation  due  solely  to  change  in  hydraulic 

boundary  conditions. 


Hydraulic  boundary  condition 
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(a)  Levee  without  impervious  riverside  membrane. 


Au 


Pore  pressure  calculated  from 
steady  state  seepage  analysis 


Time 


► 


(b)  Change  in  pore  pressure  due  to  hydraulic  boundary  conditions. 


The  true  pore  pressure  acting  at  Point  A  for  a  real  levee  at  any  given  time 
would  be  a  combination  of  that  shown  in  Figures  8.1  and  8.2.  If  the  levee 
was  constructed  of  clay,  the  pore  pressure  increase  from  the  change  in 
boundary  stress  would  dominate  for  the  initial  time  periods.  At  later  time 
periods,  the  pore  pressure  changes  from  seepage  would  be  important.  If 
the  levee  was  constructed  of  sand,  the  pore  pressure  developed  from  the 
change  in  boundary  stress  would  be  negligible,  and  the  change  in  pore 
pressure  would  mainly  be  due  to  seepage. 

In  order  for  an  effective  stress  stability  analysis  to  be  correct,  the  pore 
pressures  on  the  failure  plane  must  be  accurately  determined. 
Accommodating  both  mechanisms  of  pore  pressure  development  requires 
a  coupled  analysis.  Although  there  have  been  research  publications 
(Berilgen  2007;  Pinyol  et  al.  2008)  reporting  success  with  coupled 
analyses  for  some  types  of  problems,  these  analyses  are  not  part  of  current 
geotechnical  engineering  practice. 
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9  General  Recommendations 

Prior  to  undertaking  a  transient  seepage  analysis,  the  engineer  should 
realize  that  only  approximate  results  will  be  obtained.  Even  if  all 
parameters  necessary  to  characterize  the  soils  are  tediously  measured  in  a 
high-quality  geotechnical  laboratory,  the  results  of  the  analysis  still  may 
deviate  from  reality  (Sleep  2011). 

The  engineer  should  recognize  that  transient  seepage  analyses  only 
calculate  the  changes  in  pore  pressure  due  to  changes  in  the  hydraulic 
boundary  conditions.  The  analyses  do  not  calculate  the  changes  in  pore 
pressures  due  to  changes  in  total  stress  or  shear  stress.  The  actual  field 
pore  pressure  may  differ  considerably  from  that  calculated  using  a 
transient  seepage  analysis. 

General  recommendations  are  as  follows: 

1.  Extreme  caution  should  be  used  if  the  results  of  the  transient  seepage 
analysis  are  to  be  used  in  a  slope  stability  analysis.  The  pore  pressures 
calculated  will  not  be  correct  unless  all  pore  pressures  generated  from  the 
boundary  normal  and  shear  stresses  are  dissipated. 

2.  If  the  engineer  does  not  have  experience  in  using  the  particular  seepage 
program,  the  i-D  example  problems  and  the  generic  levee  section  included 
in  this  report  should  be  analyzed  first.  Only  after  equivalent  results  are 
obtained  should  the  user  attempt  the  program  on  an  engineering  project. 

3.  Although  all  of  the  soil  property  values  are  important,  the  user  should 
understand  that  the  compressibility  has  a  significant  effect  on  the  changes 
in  pore  pressure  over  time.  If  the  compressibility  is  equal  to  zero,  which  is 
sometimes  the  default  value  for  some  computer  programs,  the  transient 
analysis  is  essentially  the  equivalent  of  a  series  of  steady-state  analyses  for 
different  values  of  the  total  head  from  the  hydrograph  if  the  flow  is  only 
saturated.  The  greater  the  value  of  compressibility,  the  slower  the  pore 
pressures  change  relative  to  the  change  in  hydraulic  boundary  conditions. 

4.  Special  attention  should  be  paid  to  the  convergence  of  the  finite  element 
solution.  Users  should  acquaint  themselves  with  the  features  in  the 
particular  program  that  allow  the  user  to  determine  if  the  solution 
converged.  This  should  be  done  on  every  run.  An  unconverged  solution  is 
not  reliable.  This  is  especially  critical  in  transient  computations  because 
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errors  at  a  given  time  will  propagate  to  the  next  time  step  and,  in  fact,  all 
time  steps  to  follow. 

5.  The  use  of  two  different  computer  programs  may  be  helpful  in  identifying 
errors  in  the  analysis.  If  the  analyses  are  conducted  using  only  one 
program,  many  errors  may  go  undetected. 

6.  The  hydraulic  boundary  conditions  used  in  a  transient  analysis  are  just  as 
important  as  they  are  in  a  steady-state  analysis.  Recommendations 
regarding  selecting  the  appropriate  hydraulic  boundary  conditions  can  be 
found  in  Brandon  et  al.  (2014). 
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Appendix  A;  Derivation  of  the  1-D 
Unsaturated  Fiow  Soiution 

The  derivation  of  the  i-D  unsaturated  flow  problem  presented  in 
Chapter  6  is  patterned  after  Tracy  (2011)  and  described  in  the  following 
sections. 

A.l  Description  of  the  probiem 

As  described  in  Chapter  6  and  illustrated  in  Figure  6.1,  a  column  of  soil  of 
height,  L,  is  initially  almost  dry  with  a  total  head  oiht  —  Pr  +  z  or 
equivalently  a  pressure  head  of  p  =  p,.  until  a  pool  of  water  from 
infiltration  is  applied  at  the  top  (z  =  L)  such  that  the  total  head  remains 
ht  =  L  for  a  period  of  time  at  z  =  L.  The  boundary  condition  at  the  bottom 
(z  =  0)  remains  hf-  =  Pr- 

A.2  Hydrauiic  conductivity 

Hydraulic  conductivity  (/c)  is  represented  by  a  saturated  hydraulic  conduc¬ 
tivity  (kg)  and  a  multiplication  factor,  0  <  kr  <  1,  often  referred  to  as  rela¬ 
tive  hydraulic  conductivity.  Thus, 


k 


KK 


(A.l) 


A.2.1  Gardner’s  equation 

Finding  analytical  solutions  for  unsaturated  flow  is  very  difficult.  One 
approach  is  to  use  a  simpler  function  for  relative  hydraulic  conductivity 
than,  for  instance,  the  van  Genuchten  equation  (1980).  Gardner’s  equation 
(1958)  is  used  and  given  below. 

K  =  e^\p<0  (A.2) 


where: 


a  =  positive  parameter  (i/L) 
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A.3  Moisture  content 

Volumetric  moisture  content,  9,  as  computed  by  van  Genuchten  (1980) 
equation  and  shown  in  Figure  7.14,  is  too  complicated  to  be  used  in  this 
analytic  solution  derivation,  so  a  simpler  expression  is  used.  In  fact,  what 
works  is  that  dk/dO  be  a  constant  (Warrick  2003),  which  converts  Equa¬ 
tion  3.5  to  a  linear  PDE  when  using  Equation  A.2.  This  constant  derivative 
requirement  is  accomplished  here  by  first  representing  9  in  terms  of  the 
effective  saturation,  Sg. 


S  = 


O-Or 

9-9. 


(A.3) 


where: 

9r  =  residual  moisture  content 
9s  =  saturated  moisture  content 
Sg  =  effective  saturation 

Next  since  0  <  5g  <  1  just  like  the  simpler  ky  function,  Sg  is  approximated 


(A.4) 


The  moisture  content  becomes 


9  =  K{9,-9,)  +  9, 


(A.5) 


A.4  Governing  partial  differential  equation  (PDE) 

Equation  3.5  reduces  to  the  following  equation  for  i-D  unsaturated  flow  in 
a  homogeneous,  isotropic  medium  with  no  sources  or  sinks: 


A 

dz 


dt 


(A.6) 


It  is  more  convenient  to  use  a  PDE  in  terms  of  pressure  head  (p),  so 
Equation  A.6  becomes 
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A 

dz 


kk 


d{p  +  z) 
dz 


dt 


JL 

dz 


A 

dz 


dp 

dz 


+  1 


K  dt 


ijc  — ]  I  ^  ^{K{^s~^r)  +  ^r 

’’  dz  i 


dz 


dt 


d 

\k  ^P] 

1  PK  _ 

dk^ 

dz 

roz] 

dz 

[  K  J 

dt 

(A.7) 


A.5  First  change  of  variables 

Equation  A.7  is  a  nonlinear  PDE.  The  first  challenge  is  to  convert  this 
equation  to  a  linear  PDE.  This  is  done  by  the  change  of  variables, 

p  =  e“^-e  (A.8) 


where: 


e  = 


e 


Pr  =  pressure  head  when  the  soil  is  at  residual  moisture  conditions 
Using  Equations  A.2,  A.7,  and  A.8  gives 

^  =  (A.c 

dz  a  dz 


A 

dz 


dp 


1 

e""'— e 

a  dz 


d[p  +  e 


dz 


(^s-^r)^(P  +  g) 

K  9t 


1  d^p  ^  dp  _  {9^-9,)  dp 
a  dz^  dz  k^  dt 


d^p  dp  dp 
dz^  dz  dt 


(A.IO) 


where: 
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C 


a(9.-e.) 

K 


Equation  A.io  is  the  new  PE  using  the  variable,  p. 

A.6  Initial  conditions 

The  initial  conditions  for  Equation  A.io  are 

p(z,0)  =  —  e  =  0 

A.7  Boundary  conditions 

The  boundary  conditions  for  Equation  A.io  are 

p(0,f)  =  0 


p{L,t)  =  - e  =  1  - e 


A.8  Steady-state  solution 

The  governing  PDE  for  the  steady-state  solution  is 


9^Pss  I  ^dp^s 
dz^  dz 


0 


The  general  solution  is 


=  a^e 


a^e 


where: 


%  =  constant  evaluated  from  boundary  conditions 
a2  =  constant  evaluated  from  boundary  conditions 

with  rui  and  m2  being  solutions  to 

+am  =  0 


(A.ll) 


(A.12) 

(A.13) 


(A.  14) 


(A.15) 


(A.16) 


The  two  solutions  to  Equation  A.16  are 
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m  =  0,-a 


and  Equation  A.  15  becomes 


(A.17) 


From  boundary  conditions, 


0  =  +  02 

1  -  e  =  +  a^e 


-oL 


02  = 


1-e 


1-e 


-aL 


}CL±  ^2 


(A.18) 


An  alternate  form  that  will  be  used  here  is 


—z  — z 
>2  _/:>  2 


-e  2 


-  /  \  -(i- 

Pss=(l-£)e' 


sinh 


a 


sinh 


2 


(A.19) 


A.9  Second  change  of  variables 

To  solve  Equation  A.10,  a  second  change  of  variables  is  used.  It  is 


P  =  P-Ps 


(A.20) 


The  new  PDF  starting  with  Equation  A.10  becomes 

d^{P  +  Pss)  ^  ^^(P  +  PsJ_^^(P  +  Pss) 


dz 


dz 


dt 
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a- 


9Pss' 

1  1 

dz 

) 

dz^ 

dz  dt 


r\0  ^ 

^.  +  a^  =  c^ 

dz^  dz  dt 


(A.21) 


A.IO  Initial  conditions 

The  initial  conditions  for  Equation  A.21  using  Equation  A.11  and  A.20  are 


p(z,0)  =  -p, 


(A.22) 


A.ll  Boundary  conditions 

Because  the  boundary  conditions  for  Equation  21  are  the  same  as  the 
steady-state  boundary  conditions,  the  boundary  conditions  for 
Equation  A.21  using  Equation  A.20  become 

p(0,f)  =  0  (A.23) 


p{L,t)  =  0 


(A.24) 


A.12  Separation  of  variables 

A  common  way  to  solve  Equation  A.21  is  to  employ  the  use  of  separation  of 
variables.  This  is  done  by  assuming  a  solution  of  the  form, 

p(z,f)  =  ^(z)T(f)  (A.25) 


Equation  A.25  is  substituted  into  Equation  A.21  to  produce 


“I”  ^ 


dz^ 


dz 


dt 


.  .5^[^(z)]  5[^(z)]  dlrit 


dz 


dz 


dt 


(A.26) 


Shortening  the  notation  and  realizing  that  the  partial  derivatives  can  now 
be  converted  to  total  derivatives.  Equation  A.26  becomes 
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dr 


^  far— =  ct  — 
dz  dz  dt 


Dividing  both  sides  of  Equation  A.  27  by  gives 


1 

\d^^  d^l 

,  +a 

1 

=  c— 

dr 

dz^  dz 

\  J 

T 

dt, 

(A.27) 


(A.28) 


Each  side  of  the  =  sign  must  be  equal  to  some  constant,  say  rj,  to  allow  for 
a  solution  to  exist.  This  produces  the  two  equations  that  must  be  solved 
independently. 


d^ 

dz^  dz 


=  rj,c- 


dr 


=  1 


(A.29) 


Eor  convenience,  choose 


ttIc 

^.=-^-A^A,=^,A:  =  0,1,2,. 


(A.30) 


This  gives 


d^^  d^ 
dz^  ^“dz^ 


r  r^dT  1 

^  =  0— +  - 

(  2  'i 

4  I 

\  J 

dt  c 

1  4  J 

T  =  0 


(A.31) 


Using  the  same  procedure  as  before. 


^=0  2  0,0"^^"+ V 


,i  =  V^,T  =  0 


(A.32) 


where: 


4 


Ufc  =  constant  evaluated  from  boundary  conditions 
=  constant  evaluated  from  boundary  conditions 


It  is  more  convenient  to  use 
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A  sin(A,z)  + 5,  COS  (A,z) 


(A.33) 


where: 

i4fc  =  constant  evaluated  from  boundary  conditions 
B,^  =  constant  evaluated  from  boundary  conditions 

Boundary  conditions  require  5/^  =  0.  The  solution  for  p  now  takes  the 
form, 


P  =  LAsin(A,z)e 


k=l 


(A.34) 


From  initial  conditions, 


sinh 

a 

—z 

[2 

sinh 

'^L 

,2  , 

e  2  ^Aj.sin(A^z)  (A.35) 


fc=i 


Using  Fourier  series  evaluation. 


sinh 


a 


L  0 


^  /  A 

/Of 

sinh  — z  sin(Ajtz)(iz  (A. 36) 


The  integral  can  be  evaluated  by  integration  by  parts  to  obtain 

L  /  \ 

1  =  J'  sinh  — z  sin(Aj.z)dz 


=  -1 


k+l  A. 


-sinh 


+  A^. 


2 


I  =  (— l)''^^^^sinh 
CBk 


(A.37) 
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The  solution  for  p  is  therefore 

p  ^  V (- 1)"^  f ^1  sin  ( A,z)  (A. 38) 

Lc  tt"  '  [Hkj 

From  Equations  A.8,  A.19,  A.20,  A,  30,  and  A.38,  the  final  solution  is 


P  =  P  +  Pss 

(A.39) 

p  =  -ln(p  +  e) 

(A.40) 

a 

h^=p  +  z 

(A.41) 
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Appendix  B;  Derivation  of  the  1-D  Saturated 
Fiow  Soiution 

This  appendix  gives  a  derivation  for  the  i-D  saturated  flow  problem  given 
in  Chapter  6  and  is  patterned  after  the  derivation  in  Appendix  A. 

B.l  Description  of  the  probiem 

A  totally  saturated  column  of  soil  of  height,  L,  as  shown  in  Figure  6.1  is 
initially  at  hydrostatic  conditions  of  h^-  =  Hq  =  L  until  a  total  head  of 
hf-  =  is  applied  at  the  top  {z  -  L)  and  remains  at  this  level  for  an 
indefinite  period  of  time.  The  boundary  condition  at  the  bottom  (z  =  0) 
remains  hf-  =  Hq. 

B.2  Governing  PDE 

Equation  3.5  reduces  to  the  following  equation  for  a  saturated,  homoge¬ 
neous,  isotropic  soil: 


c)^hj  _  1  f}h^  _ 
dz^  dt  ’  Y^m^ 

where: 

Cy  =  coefficient  of  consolidation  (L2/T) 

B.3  Steady-state  soiution 

The  steady-state  version  of  Equation  B.i  is  simply 


The  general  solution  to  Equation  B.2  is 


(B.l) 


(B.2) 


htss  =  ai+a2Z 


(B.3) 
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where: 

%  =  constant  evaluated  from  boundary  conditions 
a2  =  constant  evaluated  from  boundary  conditions 

From  boundary  conditions, 

+  ^°z  (B.4) 


B.4  Change  of  variables 

To  solve  Equation  B.i,  a  change  of  variables  is  used  as  follows: 

h,  =  h,-h,..  (B.5) 

The  new  PDE  becomes 


La 

^^hj  +h, 

dz^ 

Cv 

dt 

d^Ks  , 

a'hj 

_  1  dh, 

dz^ 

dz 

dt 

Using  Equation  B.2  with  the  understanding  that  h^ss  is  a  function  only  of  z 
gives, 


5^hj  _  1  5h^ 
dz^  fh: 

The  initial  conditions  for  Equation  B.6  are 


hj(z,0)  =  Ho-hj^^ 


(B.6) 


(B.7) 


and  the  boundary  conditions  are 


ht(0,t)  =  0 


(B.8) 


ERDC  TR-16-8 


97 


ht(L,t)  =  0 


B.5  Separation  of  variables 


A  solution  of  Equation  B.6  is  assumed  to  take  the  form, 


h^(z,t)  =  ?(z)T(t 


(B.IO) 


Putting  Equation  B.io  into  Equation  B.6  produces 


^(z)T(t)  1  d  ?(z)T(t) 


c„  dt 


(B.ll) 


dz^  c,.  dt 


(B.12) 


Dividing  by  gives 


l[d^^]_  1  [dx^ 

ddzH  c^x[dt, 


(B.13) 


Each  side  of  the  =  sign  must  be  equal  to  a  constant,  say  q,  which  yields 


?  dz^ 


1  [dx 
c„x  dt 


^  TiT  =h 


(B.14) 


Eor  convenience,  choose 


n  =  -x^,  x,  =  ^,k  =  o,i,2,. 


(B.15) 


Equation  B.14  now  becomes 


The  individual  solutions  are 
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5  =  A^sin(X^z)  +  B^cos(X^z),  T  =  e“''''*  (B.16) 

where: 

i4fc  =  constant  evaluated  from  boundary  conditions 
=  constant  evaluated  from  boundary  conditions 

The  boundary  condition  at  z  =  0  requires  =  0.  The  solution  for  ht  now 
takes  the  form, 

h^=^A,sin(A,z)e-"‘‘^  (B.17) 

k=l 

From  initial  conditions, 

TT  _  TT  (X) 

- — ^z  =  ^AkSin(Akz)  (B.18) 

-L'  k=l 
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B.5 


Final  solution 

The  final  solution  is  obtained  by  first  solving  for  from  Equation  B.5  as 
follows: 


hj  =h^  +h^ 


(B.21) 


Next,  Equation  B.4  for  and  Equation  B.20  for  ht  are  substituted  into 
Equation  B.21  to  obtain  the  final  solution, 


h,=H„  + 


H,-  H„ 


k=l 


sin(X]^z)e^ 


(B.22) 
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